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Sum rules for Jacobi matrices 
and their applications to spectral theory 

By Rowan Killip and Barry Simon* 


Abstract 

We discuss the proof of and systematic application of Case’s sum rules 
for Jacobi matrices. Of special interest is a linear combination of two of his 
sum rules which has strictly positive terms. Among our results are a complete 
classification of the spectral measures of all Jacobi matrices J for which J — Jo 
is Hilbert-Schmidt, and a proof of Nevai’s conjecture that the Szegd condition 
holds if J — Jo is trace class. 


1. Introduction 


In this paper, we will look at the spectral theory of Jacobi matrices, that 
is, infinite tridiagonal matrices, 
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with ttj > 0 and bj G R. We suppose that the entries of J are bounded, that is, 
sup„ I On I +sup„ | 6 ri| < oo so that J defines a bounded self-adjoint operator on 
= ^^({1,2,...}). Let 5j be the obvious vector in that is, with 

components 5jn which are 1 if n = j and 0 if n 7 ^ j. 

The spectral measure we associate to J is the one given by the spectral 
theorem for the vector Ji. That is, the measure /U defined by 


( 1 . 2 ) 
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There is a one-to-one correspondence between bounded Jacobi matrices 
and unit measures whose support is both compact and contains an infinite 
number of points. As we have described, one goes from J to /x by the spectral 
theorem. One way to find J, given /x, is via orthogonal polynomials. Apply¬ 
ing the Gram-Schmidt process to one gets orthonormal polynomials 

Pn{x) = + ■ ■ ■ with > 0 and 

(1.3) J Pni^x)Pyni^X^ — dfrirri' 

These polynomials obey a three-term recurrence: 

(1.4) xPyi(x^ — (ln-\-iPn-{-l{x') J- bji+lPnix') J- CLnPn—lix'), 


where Un, bn are the Jacobi matrix coefficients of the Jacobi matrix with spec¬ 
tral measure fj, (and P_i = 0 ). 

The more usual convention in the orthogonal polynomial literature is to 
start numbering of {an} and {bn} with n = 0 and then to have (1.4) with 
{an,bn,an-i) instead of (on+i, 6 ^+ 1 , a^). We made our choice to start num¬ 
bering of J at n = 1 so that we could have z"' for the free Jost function (well 
known in the physics literature with z = e*^) and yet arrange for the Jost 
function to be regular at z = 0. (Case’s Jost function in [ 6 , 7] has a pole since 
where we use uq below, he uses tt_i because his numbering starts at n = 0 .) 
There is, in any event, a notational conundrum which we solved in a way that 
we hope will not offend too many. 

An alternate way of recovering J from fj, is the continued fraction expan¬ 
sion for the function m^{z) near infinity. 


(1.5) 




-E + hi 


1 


-E + h2 + --- 


Both methods for finding J essentially go back to Stieltjes’ monumental 
paper [57]. Three-term recurrence relations appeared earlier in the work of 
Chebyshev and Markov but, of course, Stieltjes was the first to consider general 
measures in this context. While [57] does not have the continued fraction 
expansion given in (1.5), Stieltjes did discuss (1.5) elsewhere. Wall [62] calls 

(1.5) a J-fraction and the fractions used in [57], he calls S'-fractions. This has 
been discussed in many places, for example, [24], [56]. 

That every J corresponds to a spectral measure is known in the orthog¬ 
onal polynomial literature as Favard’s theorem (after Favard [15]). As noted, 
it is a consequence for bounded J of Hilbert’s spectral theorem for bounded 
operators. This appears already in the Hellinger-Toeplitz encyclopedic arti¬ 
cle [26]. Even for the general unbounded case, Stone’s book [58] noted this 
consequence before Favard’s work. 
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Given the one-to-one correspondence between /x’s and J’s, it is natural to 
ask how properties of one are reflected in the other. One is especially interested 
in J’s “close” to the free matrix, Jq with = 1 and bn = 0, that is, 


( 1 . 6 ) 


-^0 


/ 0 1 0 0 ... \ 

1 0 1 0 ... 

0 1 0 1 ... 

V 0 0 1 0 ... / 


In the orthogonal polynomial literature, the free Jacobi matrix is taken 
as ^ of our Jq since then the associated orthogonal polynomials are precisely 
Chebyshev polynomials (of the second kind). As a result, the spectral measure 
of our Jq is supported by [—2, 2] and the natural parametrization \s E = 2 cos 6. 
Here is one of our main results: 


Theorem 1. Let J he a Jacobi matrix and fi the corresponding speetral 
measure. The operator J — Jq is Hilbert-Schmidt, that is, 

(1.7) 2Y,{an-lf + Y.bl<oo 

n 

if and only if p, has the following four properties: 

(0) (Blumenthal-Weyl Criterion) The support of p, is [—2, 2] U {E^}f=i ^ 
{E~}f=i where N± are each zero, finite, or infinite, and E^ > E 2 > 

■ ■ ■ >2 and Ef < Ef < ■ ■ ■ < —2 and if N± is infinite, then 
lim^^oo E^ = ±2. 

(1) (Quasi-Szego Condition) Let fiac{E) = f{E) dE where is the Lebesgue 
absolutely continuous component of fi. Then 


r2 _ 

(1.8) J log[f{E)] \/4 -E^dE> - 00 . 

(2) (Lieb-Thirring Bound) 

N+ N- 

(1.9) ^|E+-2|3/2 + ^|^- + 2|3/2 < 00 . 

1=1 1=1 

(3) (Normalization) / dyL{E) = 1. 

Remarks. 1. Condition (0) is just a quantitative way of writing that the 
essential spectrum of J is the same as that of Jq, viz. [—2,2], consistent with 
the compactness of J — Jq. This is, of course, Weyl’s invariance theorem [63], 
[45]. Earlier, Blumenthal [5] proved something close to this in spirit for the 
case of orthogonal polynomials. 

2. Equation (1.9) is a Jacobi analog of a celebrated bound of Lieb and 
Thirring [37], [38] for Schrodinger operators. That it holds if J — Jq is Hilbert- 
Schmidt has also been recently proven by Hundertmark-Simon [27], although 
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we do not use the |-bound of [27] below. We essentially reprove (1.9) if (1.7) 
holds. 

3. We call (1.8) the quasi-Szegd condition to distinguish it from the Szego 
condition, 

r2 

( 1 . 10 ) j log[f{E)]{A-E‘^)-^/‘^dE>-oo. 


This is stronger than (1.8) although the difference only matters if / vanishes 
extremely rapidly at ±2. For example, like exp(—(2 — |£'|)““) with ^ < a < |. 
Such behavior actually occurs for certain Pollaczek polynomials [ 8 ]. 


4. It will often be useful to have a single sequence ei (J), 62 (J),... obtained 


from the numbers 




T2 


by reordering so ei(J) > 62 (J) > • • • 


0 . 


By property (1), for any J with J — Jq Hilbert-Schmidt, the essential 
support of the a.c. spectrum is [—2,2]. That is, /Uac gives positive weight to 
any subset of [—2,2] with positive measure. This follows from (1.8) because 
/ cannot vanish on any such set. This observation is the Jacobi matrix ana¬ 
logue of recent results which show that (continuous and discrete) Schrodinger 
operators with potentials V G L^, p < 2, or |I7(a;)| < (1 J- a > 1/2, 

have a.c. spectrum. (It is known that the a.c. spectrum can disappear once 
p > 2 or a < 1/2.) Research in this direction began with Kiselev [29] and cul¬ 
minated in the work of Christ-Kiselev [11], Remling [47], Deift-Killip [13], and 
Killip [28]. Especially relevant here is the work of Deift-Killip who used sum 
rules for finite range perturbations to obtain an a priori estimate. Our work 
differs from theirs (and the follow-up papers of Molchanov-Novitskii-Vainberg 
[40] and Laptev-Naboko-Safronov [36]) in two critical ways: we deal with the 
half-line sum rules so the eigenvalues are the ones for the problem of interest 
and we show that the sum rules still hold in the limit. These developments are 
particularly important for the converse direction (i.e., if p obeys (0~3) then 
J — Jo is Hilbert-Schmidt). 

In Theorem 1, the only restriction on the singular part of p on [—2,2] 
is in terms of its total mass. Given any singular measure //sing supported on 
[—2, 2] with total mass less than one, there is a Jacobi matrix J obeying (1.7) 
for which this is the singular part of the spectral measure. In particular, there 
exist Jacobi matrices J with J — Jq Hilbert-Schmidt for which [—2,2] simul¬ 
taneously supports dense point spectrum, dense singular continuous spectrum 
and absolutely continuous spectrum. Similarly, the only restriction on the 
norming constants, that is, the values of p{{E^}), is that their sum must be 
less than one. 

In the related setting of Schrodinger operators on M, Denisov [14] has 
constructed an potential which gives rise to embedded singular continuous 
spectrum. In this vein see also Kiselev [30]. We realized that the key to 
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Denisov’s result was a sum rule, not the particular method he used to construct 
his potentials. We decided to focus first on the discrete case where one avoids 
certain technicalities, but are turning to the continuum case. 

While (1.8) is the natural condition when J — Jq is Hilbert-Schmidt, we 
have a one-directional result for the Szego condition. We prove the following 
conjecture of Nevai [43]: 

Theorem 2. If J — Jq is in trace class, that is, 

(1.11) ^ jOn - 1| ^ | 6 n| < OO, 

n n 

then the Szego condition (1.10) holds. 

Remark. Nevai [42] and Geronimo-Van Assche [22] prove the Szego con¬ 
dition holds under the slightly stronger hypothesis 

^(logn) la„ - 1] + ^(logn) \hn\ < oo. 

n n 

We will also prove 

If J — Jq is compact and 

I ^ 1/2 I 1/2 

Hi®;-/ +i:|®7+2| 

i j 

(ii) limsup^^j,o oi... oat > 0 
then (1.10) holds. 

We will prove Theorem 2 from Theorem 3 by using a ^ power Lieb-Thirring 
inequality, as proven by Hundertmark-Simon [27]. 

For the special case where // has no mass outside [—2,2] (i.e., = A"_ 

= 0), there are over seventy years of results related to Theorem 1 with im¬ 
portant contributions by Szego [59], [60], Shohat [49], Geronomius [23], Krein 
[33], and Kolmogorov [32]. Their results are summarized by Nevai [43] as: 

Theorem 4 (Previously Known). Suppose /j, is a probability measure 
supported on [— 2 , 2 ]. The Szego condition (1.10) holds if and only if 

(i) J — Jq is Hilbert-Schmidt. 

(ii) X](®n “ 1 ) 'ffbn are {conditionally) convergent. 

Of course, the major difference between this result and Theorem 1 is 
that we can handle bound states (i.e., eigenvalues outside [— 2 , 2 ]) and the 
methods of Szego, Shohat, and Geronimus seem unable to. Indeed, as we 


Theorem 3. 

(i) 

( 1 . 12 ) 
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will see below, the condition of no eigenvalues is very restrictive. A second 
issue is that we focus on the previously unstudied (or lightly studied; e.g., it 
is mentioned in [39]) condition which we have called the quasi-Szego condition 
(1.8), which is strictly weaker than the Szego condition (1.10). Third, related 
to the first point, we do not have any requirement for conditional convergence 
of En=l(«n-1) or En=lbn. 

The Szego condition, though, has other uses (see Szego [60], Akhiezer [2]), 
so it is a natural object independently of the issue of studying the spectral 
condition. 

We emphasize that the assumption that fi has no pure points outside 
[—2,2] is extremely strong. Indeed, while the Szego condition plus this as¬ 
sumption implies (i) and (ii) above, to deduce the Szego condition requires 
only a very small part of (ii). We 

Theorem T. If (t{J) C [-2,2] and 
(i) limsupArEn=ilog(«n) > -oo, 

then the Szego condition holds. If cr{J) C [—2, 2] and either (i) or the Szego 
condition holds, then 

(ii) E“i(an-l)2 + E“ift^<oo, 

(hi) limTv^oo log(an) exists {and is finite), 

(iv) \im]\f^ooJ2n=ibn exists {and is finite). 

In partieular, if a{ J) C [—2,2], then (i) implies (ii)-(iv). 

In Nevai [41], it is stated and proven (see pg. 124) that jon — 1| < oo 
implies the Szego condition, but it turns out that his method of proof only 
requires our condition (i). Nevai informs us that he believes his result was 
probably known to Geronimus. 

The key to our proofs is a family of sum rules stated by Case in [7]. Case 
was motivated by Flaschka’s calculation of the first integrals for the Toda 
lattice for finite [16] and doubly infinite Jacobi matrices [17]. Case’s method 
of proof is partly patterned after that of Flaschka in [17]. 

To state these rules, it is natural to change variables from E to z via 

(1.13) E = z + ^. 

We choose the solution of (1-13) with \z\ < 1, namely 

e-Ve^-a], 

where we take the branch of with > 0 for p > 0. In this way, E z is 


(1.14) 
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the conformal map of { 00 } U C\[— 2, 2]io D = {z\ | 2 ;| < 1}, which takes 00 to 
0 and (in the limit) ±2 to ±1. The points E G [—2, 2] are mapped to 2 : = 
where E = 2 cos 9. 

The conformal map suggests replacing by 

(1.15) Mf,{z) = -mf,(^E{z)'^ = -m^(^z + = J ^ ' 

We have introduced a minus sign so that ImM^( 2 ;) > 0 when Imz > 0. Note 
that ImE > 0 ^ m^(£') > 0 but E ^ z maps the upper half-plane to the 
lower half-disk. 

If jj, obeys the Blumenthal-Weyl criterion, is meromorphic on D with 
poles at the points where 

(1.16) | 7 j| > 1 and Ef = -ff + ( 7 ^)“^- 

As with E^, we renumber to a single sequence \j3i\ > |/ 32 | > • • • > 1. 

By general principles, has boundary values almost everywhere on the 
circle, 

(1.17) = hm M^{rE^) 

rtl 


with M^(e *®) = M^(e*®) and ImM^(e*^) > 0 for 6* G (0,vr) 
From the integral representation (1.2), 


(1.18) 


Imm^(Fi -|- iO) = TT 


dE 


so using dE = —2sm9d6 = —(4 — d9, the quasi-Szego condition (1.8) 

becomes 

4 [ log[ImM„(e*®)] sin^ > —00 

Jo 

and the Szego condition (1.10) is 

[ log[ImM„(e*®)] > — 00 . 

Jo 

Moreover, we have by (1.18) that 

(1.19) f [\m[M^{E^)]sm9de = n^c{-2,2)<l. 

With these notational preliminaries out of the way, we can state Case’s 
sum rules. For future reference, we give them names: 


Co: 


( 1 . 20 ) 


1 r 
dvr 7_^ 


log 


sin0 


ImM(e*®)J 


d9 = J2^og\Pj\ - ^ log I 
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and for n = 1,2,.. 


C„: 


( 1 . 21 ) 


1 

271 J-T, 


-IT- / log 


sin0 


[lmM(e*®) J 


cos{n6) dO + - /?• "■) 




where is the Chebyshev polynomial (of the first kind). 

We note that Case did not have the compact form of the right side of 
(1.21), but he used implicitly defined polynomials which he did not recognize 
as Chebyshev polynomials (though he did give explicit formulae for small n). 
Moreover, his arguments are formal. In an earlier paper, he indicates that the 
conditions he needs are 


( 1 . 22 ) \an-l\ + \bn\ <C{l + n‘^)-^ 

but he also claims this implies N+ < oo, N- < oo, and, as Chihara [9] noted, 
this is false. We believe that Case’s implicit methods could be made to work 
if — 1-1 + \bn\] < oo rather than (1.22). In any event, we will provide 

explicit proofs of the sum rules—indeed, from two points of view. 

One of our primary observations is the power of a certain combination of 
the Case sum rules, Cq + ^C 2 - It says 

P2: 


j j 

where G{a) = a? — 1 — log \af and F{E) = | [f3^ —— log ], with f3 given 
hy E = P + \P\ > 1 (cf. (1.16)). 

As with the other sum rules, the terms on the left-hand side are purely 
spectral—they can be easily found from those on the right depend in a 
simple way on the coefficients of J. 

The significance of (1-23) lies in the fact that each of its terms is non¬ 
negative. It is not difficult to see (see the end of §3) that E{E) > 0 for 
E £ W \ [—2,2] and that G{a) > 0 for a G (0, 00 ). To see that the integral is 
also nonnegative, we employ Jensen’s inequality. Notice that y ^ — log{y) is 
convex and ^ Jq sin^ 0d0 = 1 so 
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sin(0) 


sin^ 9 (16 


1 2 
2 T 



-log 


'Im M' 
sin 6 


sin^ (9) (16 


> —log — [ {Im M) sm{9) d9 
2 [tt Jo 

= -^log[^ac(-2,2)] > 0 

by (1.19). 

The hard work in this paper will be to extend the sum rule to equalities 
or inequalities in fairly general settings. Indeed, we will prove the following; 


Theorem 5. If J is a Jacobi matrix for which the right-hand side of 
(1.23) is finite, then the left-hand side is also finite and LHS < RHS. 

Theorem 6. If fi is a probability measure that obeys the Blumenthal- 
Weyl criterion and the left-hand side of (1.23) is finite, then the right-hand 
side of (1.23) is also finite and LHS > RHS. 


In other words, the P 2 sum rule always holds although both sides may 
be infinite. We will see (Proposition 3.4) that G{a) has a zero only at a = 1 
where G{a) = 2{a — 1)^ + 0{{a — 1)^) so the RHS of (1.23) is finite if and only 
if Z) — 1)^ < 00 , that is, J is Hilbert-Schmidt. On the other hand, 

we will see (see Proposition 3.5) that F{Ej) = {\Ej\ — 2 )^^ -|- 0{{\Ej\ — 2)^) 
so the LHS of (1.23) is finite if and only if the quasi-Szego condition (1.8) and 
Lieb-Thirring bound (1.9) hold. Thus, Theorems 5 and 6 imply Theorem 1. 

The major tool in proving the Case sum rules is a function that arises in 
essentially four distinct guises: 


(1) The perturbation determinant defined as 


(1.25) 


L{z‘, J) = det 


{J — z — z ^){Jo — z — z ^) ^ 


(2) The dost function, uo{z; J) defined for suitable z and J. The dost solution 
is the unique solution of 

(1.26) anUn+l + hnUn + On-lUn-l = {z + Z~^)Un 

n > 1 with ao = 1 which obeys 

(1.27) lim z~^Un = 1. 

n^oo 

The dost function is uo{z; J) = uq. 

(3) Ratio asymptotics of the orthogonal polynomials Pn, 


(1.28) 


lim Pn{z + z-^)z^. 
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(4) The Szego function, normally only defined when A^_|_ = N- = 0: 


(1.29) D{z) = exp^— J log 27rsin(0)/(2cos0) 


+ z 
— z 


d9 


where dfj, = f{E)dE + d/Usmg- 

These functions are not all equal, but they are closely related. L(z; J) 
is defined for \z\ < 1 by the trace class theory of determinants [25], [53] so 
long as J — Jo is trace class. We will see in that case it has a continuation to 
{z I |z| < 1, z 7 ^ ±1} and, when J — Jq is finite rank, it is a polynomial. The 
dost function is related to L by 


(1.30) 


no(z; J) 



L(z;J). 


Indeed, we will define all Un by formulae analogous to (1.30) and show that they 
obey (1.26)/(1.27). The dost solution is normally constructed using existence 
theory for the difference equation (1.26). We show directly that the limit in 
(1.28) is tto(J, z)/(l — z^). Finally, the connection of D{z) to uo{z) is 

(1.31) D{z) = ( 2 )- 1/2 _ ^ 2 ) jyl_ 


Connected to this formula, we will prove that 


(1.32) 



sin0 

ImM^(0) ’ 


from which (1.31) will follow easily when J — Jq is nice enough. The result for 
general trace class J — Jq is obviously new since it requires Nevai’s conjecture 
to even define D in that generality. It will require the analytic tools of this 
paper. 

In going from the formal sum rules to our general results like Theorems 4 
and 5, we will use three technical tools; 


(1) That the map ^ ^ log( ) sin^ 6 dO and the similar map with 
sin^ 6 dO replaced by dO is weakly lower semicontinuous. As we will see, 
these maps are essentially the negatives of entropies and this will be a 
known upper semicontinuity of an entropy. 

(2) Rather than prove the sum rules in one step, we will have a way to prove 
them one site at a time, which yields inequalities that go in the opposite 
direction from the semicontinuity in (1). 

(3) A detailed analysis of how eigenvalues change as a truncation is removed. 
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In Section 2, we discuss the construction and properties of the pertur¬ 
bation determinant and the dost function. In Section 3, we give a proof of 
the Case sum rules for nice enough J — Jq in the spirit of Flaschka’s [16] and 
Case’s [7] papers, and in Section 4, a second proof implementing tool (2) above. 
Section 5 discusses the Szego and quasi-Szegd integrals as entropies and the 
associated semicontinuity, and Section 6 implements tool (3). Theorem 5 is 
proven in Section 7, and Theorem 6 in Section 8. 

Section 9 discusses the Cq sum rule and proves Nevai’s conjecture. 

The proof of Nevai’s conjecture itself will be quite simple—the Cq sum 
rule and semicontinuity of the entropy will provide an inequality that shows 
the Szego integral is finite. We will have to work quite a bit harder to show that 
the sum rule holds in this case, that is, that the inequality we get is actually 
an equality. 

In Section 10, we turn to another aspect that the sum rules expose: the 
fact that a dearth of bound states forces a.c. spectrum. For Schrodinger op¬ 
erators, there are many V’s which lead to cj(— A + V) = [0,oo). This always 
happens, for example, if V{x) >0 and lim|j,|^(,o I/(x) = 0. But for discrete 
Schrodinger operators, that is, Jacobi matrices with = 1, this phenomenon 
is not widespread because it( Jq) has two sides. Making 6n A 0 to prevent eigen¬ 
values in (-00,-2) just forces them in (2,oo)! We will prove two somewhat 
surprising results (the en{J) are defined in Remark 6 after Theorem 1). 

Theorem 7. If J is a Jacobi matrix with a„ = 1 and Yin \ ^n{J)\^^‘^ < oo, 
then fTac(J) = [-2,2]. 

Theorem 8. Let W be a two-sided Jacobi matrix with = 1 and no 
eigenvalues. Then bn = 0, that is, W = Wq, the free Jacobi matrix. 

We emphasize that Theorem 8 does not presuppose any reflectionless con¬ 
dition. 

Acknowledgments. We thank F. Gesztesy, N. Makarov, P. Nevai, 
M. B. Ruskai, and V. Totik for useful discussions. R.K. would like to thank 
T. Tombrello for the hospitality of Caltech where this work was initiated. 


2. Perturbation determinants and the Jost function 


In this section we introduce the perturbation determinant 
L{z;J) =det (^J — E{z)^ (^Jq — E{z)^ ; E{z) = z + z 


-1 


and describe its analytic properties. This leads naturally to a discussion of the 
Jost function commencing with the introduction of the Jost solution (2.63). 
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The section ends with some remarks on the asymptotics of orthogonal poly¬ 
nomials. We begin, however, with notation, the basic properties of Jq, and 
a brief review of determinants for trace class and Hilbert-Schmidt operators. 
The analysis of L begins in earnest with Theorem 2.5. 

Throughout, J represents a matrix of the form (1.1) thought of as an 
operator on t'^(Z_|_). The special case a„ = 1, 6^ = 0 is denoted by Jq and 
5J = J — Jq constitutes the perturbation. If SJ is finite rank (i.e., for large n, 
On = 1 and = 0), we say that J is finite range. 

It is natural to approximate the true perturbation by one of finite rank. 
We define Jn as the semi-infinite matrix, 

^ 6i ai 0 

62 02 

(2.1) Jn — ... bn—I On—1 

On— 1 bn 1 

1 0 1 

V 1 0 

that is, Jn has 6^ = 0 for m > n and am = 1 for m > n — 1. Notice that 
Jn — Jq has rank at most n. 

We write the nxn matrix obtained by taking the first n rows and columns 
of J (or of Jn) as Jn-F- The nxn matrix formed from Jq will be called Jo-n;F- 
A different class of associated objects will be the semi-infinite matrices 
J^^'^ obtained from J by dropping the first n rows and columns of J, that is, 

bn+l Oji+1 0 
On+l bn+2 Ofj +2 
0 0 ) 2-|_2 bn-\-3 

As the next preliminary, we need some elementary facts about Jq, the free 
Jacobi matrix. Fix z with \z\ < 1. Look for solutions of 

(2.3) Un^l J- nn —1 — (^Z -\- Z )'Un) n ‘2 

as sequences without any a priori conditions at infinity or n = 1. The solutions 
of (2.3) are linear combinations of the two “obvious” solutions given by 

(2.4) ut{z) = 

Note that is at infinity since \z\ < 1. The linear combination that obeys 

U 2 = {z + z~^)ui 





as required by the matrix ending at zero is (unique up to a constant) 
(2.5) u(J\z) = z-^ - zA 
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Noting that the Wronskian of and it"'' is z ^ — -z, we see that ( Jq — E(z)) ^ 

has the matrix elements ——^)~^^min(n,m)('^)^max(n,m)('^) either by a direct 
calculation or standard Green’s function formula. We have thus proven that 


(2.6) (Jo - E{z))^l 

(2.7) 






min(m,n) —1 

- E - 

i=o 


l+\m—n\+2j 


where the second comes from {z ^ — z)(z^ ^ + z^ " + ••• + 2 ;’^ ^ n _ 

by telescoping. (2.7) has two implications we will need later: 


( 2 . 8 ) 


Ul < 1 ^ 


{Jo — E{z)nm < niin(n, m) 


and that while the operator (Jq — E{z))~^ becomes singular as \z\ t 1, the 
matrix elements do not; indeed, they are polynomials in z. 

We need an additional fact about Jq; 


Proposition 2.1. The characteristic polynomial of Jo-n-,F is 

/ -n-l _ ^n+l\ 

(2.9) det( j;(z) - Jo,n;F) = ^ [h ^(^)) 

where Un{cos6) = sin[(n + l)0]/sin(0) is the Chebyshev polynomial of the 
second kind. In particular, 


( 2 . 10 ) 


det[.E(z) - Jo;n+i;F] _ 

™ det[E{z) - Jo;n;F] 


Proof. Let 

(2.11) gn{z) = det(j;(2;) - Jo;n;F)- 

By expanding in minors 

gn+ 2 {z) = {z + z~^)gn+i{z) - gn{z). 

Given that gi = z + z~^ and go = 1, we obtain the first equality of (2.9) by 
induction. The second equality and (2.10) then follow easily. □ 

In Section 4, we will need 

Proposition 2.2. LetTm be the Chebyshev polynomial {of the first kind): 


(2.12) 


Tm{cos6) = cos{m6). 



266 


ROWAN KILLIP AND BARRY SIMON 


Then 

n m = 2i{n + 1); £ G Z 

— ^ — ^(—I)*” otherwise. 

In particular, for m fixed, once n > — 1 the trace is independent of n. 


(2.13) 


Tr 


T 

± r 




n\F 


Proof. As noted above, the characteristic polynomial of Jo,n;F is Un{E/2). 
That is, det[2cos(0) — Jo;n;F] = sin[(n + l)0]/sin[0]. This implies that the 
eigenvalues of Jo-n-,F are given by 


(2.14) 


= 2 


cos 


/CTT 

n + 1 


A: = 1,..., n. 


So by (2.12), Tm{^E^n^) = cos (gf). Thus, 


n 

Tr r„(ijo;n;F) = XI 


cos 


k=l 


kniF \ 
n + l) 


n+l 


+ i i: -p(()^)- 

k=-n V -L / 


The final sum is 2n + 2 if m is 


a multiple of 2(n + 1) and 0 if it is not. 


□ 


As a final preliminary, we discuss Hilbert space determinants [25], [52], 
[53]. Let Ip denote the Schatten classes of operators with norm ||A||p = 
Tr(|A|^) as described for example, in [53]. In particular, Ii and I 2 are the 
trace class and Hilbert-Schmidt operators, respectively. 

For each A £ Ii, one can define a comp lex-valued function det(l -|-A) (see 
[25], [53], [52]), so that 

(2.15) |det(l + A)| < exp(||A||i) 


and A det(l -|- A) is continuous; indeed [53, pg. 48], 

(2.16) |det(l + A) - det(l -|- H)| < ||A - H||i exp(||A||i -|- ||H||i -|- 1). 
We will also use the following properties: 

(2.17) A,BgIi ^ det(l-|-A) det(l-|-H) = det(l-|-A-|-H-|-AH) 

(2.18) AB,BA£li det(l-|-AH) = det(l-|-HA) 

(2.19) (1 -|- A) is invertible if and only if det(l -|- A) / 0 

(2.20) 2 : A( 2 ;) analytic ^ det(l -|- A( 2 :)) analytic. 

If A is finite rank and H is a finite-dimensional self-adjoint projection, 

(2.21) PAP = A ^ det(l -|- A) = detp 7 ^(lp'>^ -|- PAP), 


where detp-^^ is the standard finite-dimensional determinant. 
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For A e X 2 , (1 + A)e~^ — 1 G Ji, so one defines (see [53, pp. 106-108]) 

(2.22) det 2 (l + ^) = det((l + ^)e“"^). 

Then 

(2.23) |det2(l + A)| < exp(plll) 

(2.24) |det 2 (l + A)-det 2 (l + B)| < P - 5||2 exp((P ||2 + ||5||2 + 1)^) 
and, if ^ G Ti, 

(2.25) det2(l + ^) = det(l + Al)e-^(^) 
or 

(2.26) det(l + A)= det 2 (l + 

To estimate the Zp norms of operators we use 

Lemma 2.3. If A is a matrix and || • jjp the Schatten Zp norm [53], then 

(i) 

(2.27) ll^lli = ^ la„^l^ 

n,m 

(ii) 

(2.28) Pl|i<^|an,„l, 

n,m 

(iii) For any j and p, 

(2.29) ^ |o„,< \\A\\l. 


Proof, (i) is standard, (ii) follows from the triangle inequality for jj • jji 
and the fact that a matrix which a single nonzero matrix element, a, has trace 
norm jaj. (iii) follows from a result of Simon [53], [51] that 


P = sup<j I (Pr, 


, A^'lpn) 1 ^ 


{<Pn}, {f’n} orthonormal sets 


□ 


The following factorization will often be useful. Define 
Cn = max(la„_i - 1] , J^nl , kn - 1|) 

which is the maximum matrix element in the row and column. Let C 
be the diagonal matrix with matrix elements Cn- Define U by 

(2.30) 6J = 

Then [7 is a tridiagonal matrix with matrix elements bounded by 1 so 

(2.31) PI! <3. 
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One use of (2.30) is the following: 

Theorem 2.4. Let Cn = max(|an-i — 1| , \hn\ , \an — 1|)- For any p £ 
[l,oo), 

/ \ 

(2.32) 5 < ll<J^llp<3 • 

^ n ' ^ n ^ 

Proof. The right side is immediate from (2.30) and Holder’s inequality for 
trace ideals [53]. The leftmost inequality follows from (2.29) and 

(Tki") < i;i‘'«r) +2 EK-ir ■ ^ 


With these preliminaries out of the way, we can begin discussing the per¬ 
turbation determinant L. For any J with 5J £ Zi (by (2.32) this is equivalent 
to Z) l®n — 1| + Z) \^n\ < oo), we define 


(2.33) 


L{z; J) = det [(j - E{z)'^ (jo - E{z)'^ 


-In 


for all HI < 1. Since 


(2.34) (J - E){Jo - E)-^ = 1 + 5J{Jo - E)-\ 


the determinant in (2.33) is of the form 1 + A with A £li. 

Theorem 2.5. Suppose 5J £ Ii. 

(i) L{z; J) is analytic in D = {z \ \z\ < 1}. 

(ii) L{z; J) has a zero in D only at points Zj where E{zj) is an eigenvalue 
of J, and it has zeros at all such points. All zeros are simple. 


(iii) If .J is finite range, then L{z] J) is a polynomial and so has an analytic 
continuation to all of C. 


Proof, (i) follows from (2.20). 

(ii) If Eq = E{zo) is not an eigenvalue of J, then Eq ^ cr{J) since E : 
D C\[—2,2] and aess{J) = [~2,2]. Thus, {J — Eq)/{Jo — Eq) has an inverse 
(namely, (Jq — Eq)/{J — Eq)), and so by (2.19), L{z;J) 0. If Eq is an 
eigenvalue, (J — Eq)/{Jq — Eq) is not invertible, so by (2.19), L{zq;J) = 0. 
Finally, if E{zq) is an eigenvalue, eigenvalues of J are simple by a Wronskian 
argument. That L has a simple zero under these circumstances comes from 
the following. 
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If P is the projection onto the eigenvector at Eq = E{zo), then 
(J — — P) has a removable singularity at Eq. Define 

(2.35) C{E) = {J-E)-^{1-P) + P 


so 

(2.36) {J-E)C{E) = l-P+{Eo-E)P. 

Define 


(2.37) 


D{E) = {Jo-E)C{E) 

= -6JC{E) + { J - E)C{E) 

= l-P + {Eo- E)P - SJC{E) 

= 1 + trace class. 


Moreover, 

D{E)[{J-E)/{Jo-E)] = {Jo-E)[l-P + {Eo-E)P]{Jo-E)-^ 

= l + iJo-E)[-P+iEo-E)P]iJo-E)-\ 

Thus by (2.17) first and then (2.18), 


det{D{E{z)))L{z;J) = det{l + {Jq - E)[-P + {Eq - E)P]{Jo - E)-^) 

= det(l -P + {Eo- E)P) 

= Eq- E{z), 


where we used (2.21) in the last step. Since L{z; J) has a zero at zq and 
Eq — E{z) = {z — zo)[l — has a simple zero, L{z; J) has a simple zero. 

(hi) Suppose 5J has range N, that is, N = max{n | \bn\ + |a„_i — 1| > 0} 
and let P^^') be the projection onto the span of As P^^'>5J = SJ, 

SJ{Jo - E)-^ = - E)-\ 


By (2.18), 

L{z; J) = det (l + - E{z)^ . 

Thus by (2.7), L{z; J) is a polynomial if 5J is finite range. □ 

Remarks. 1. By this argument, if 5J has range n, L{z-, J) is the determi¬ 
nant of an n X n matrix whose ij element is a polynomial of degree i + j + 1. 
That implies that we have shown L{z] J) is a polynomial of degree at most 
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2n(n + l)/2 + n = (n + 1)^. We will show later it is actually a polynomial of 
degree at most 2n — 1. 

2. The same idea shows that if Yl,n |(®n “ < oo for some 

/9 > 1, then Jo — z — is trace class for \z\ < p, and thus L{z] J) 

has an analytic continuation to {z | \z\ < p}. 

We are now interested in showing that L{z] J), defined initially only on D, 
can be continued to dD or part of dD. Our goal is to show; 

(i) If 

OO 

(2.38) ^ n[|an - 1| + |6n|] < OO, 

n=l 

then L(z; J) can be continued to all of D, that is, extends to a function 
continuous on D and analytic in D. 

(ii) For the general trace class situation, L{z-, J) has a continuation to 

(hi) As X real approaches ±1, |L(x; J)| is bounded by exp{o(l)/(l — |x|)}. 

We could interpolate between (i) and (hi) and obtain more information 
about cases where (2.38) has n replaced by n“ with 0 < a < 1 or even logn (as 
is done in [42], [22]), but using the theory of Nevanlinna functions and (hi), we 
will be able to handle the general trace class case (in Section 9), so we forgo 
these intermediate results. 

Lemma 2.6. Let C be diagonal positive trace class matrix. For \z\ < 1, 
define 

(2.39) A{z) = - E{z))-^C^/^. 

Then, as a Hilbert-Schmidt operator-valued function, A{z) extends continu¬ 
ously to D\ {—1,1}. If 

(2.40) E nCn < OO, 

n 

it has a Hilbert-Schmidt continuation to D. 

Proof. Let Anm{z) be the matrix elements of A(z). It follows from \z\ < 1 
and (2.6)/(2.8) that 

\Anmiz)\ < \Z - \Z P \ 

\Anm{z)\ < min{m,n)c]l‘^c]l^ 


(2.41) 

(2.42) 


-1 
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and each An^m{z) has a continuous extension to D. It follows from (2.41), the 
dominated convergence theorem, and 

n,m 

that so long as z stays away from ±1, {Amn{z)}n,m is continuous in the space 
£^((l,oo) X (l,oo)) so A{z) is Hilbert-Schmidt and continuous on L>\{—1,1}. 
Moreover, (2.42) and 



[min(m,n)cy^c^^ 

n,m 


< E TTlTlCfiCui 

mn 



imply that A(z) is Hilbert-Schmidt on D if (2.40) holds. 


□ 


Remark. When (2.40) holds—indeed, when 


(2.43) 


E 


n"c„ < oo 


for some a > 0—we believe that one can show A{z) has trace class boundary 
values on 5Z1\{—1,1} but we will not provide all the details since the Hilbert- 
Schmidt result suffices. To see this trace class result, we note that Imyl( 2 :) = 
{A{z) — A*{z))/2i has a rank 1 boundary value as z —> e*®; explicitly, 


(2.44) 


Imyl(e*®)mn 


—c 


1/2„V2 

n '-'m 


(sin m9){sinn6) 
{sinO) 


Thus, Imyl(e*^) is trace class and is Holder continuous in the trace norm if 
(2.43) holds. Now ReA(e*^) is the Hilbert transform of a Holder continuous 
trace class operator-valued function and so trace class. This is because when 
a function is Holder continuous, its Hilbert transform is given by a convergent 
integral, hence limit of Riemann sums. Because of potential singularities at 
±1, the details will be involved. 


Lemma 2.7. Let 6J be traee class. Then 
(2.45) t{z) = Tt{{6J){Jo - E{z))-^) 

has a continuation to l)\{—1,1}. //(2.38) holds, t{z) can be continued to D. 

Remark. We are only claiming t{z) can be continued to dD, not that 
it equals the trace of {SJ){Jo — E{z))~^ since 5J{Jo — E{z))~^ is not even a 
bounded operator for z G dD\ 

Proof. t{z) = ti{z) + t 2 {z) + t 3 {z) where 

* i (^) = J2^m-e{z))-^ 

t2{z) = ^(an-l)(Jo--E(^))El,n 

^ 3 (^) = E(«^ “ - ^(^))E+i • 
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Since, by (2.6), (2.8), 

{Jq-E{z))^^ < 2|2;-1| ^|z + l| ^ 

(Jo - £^( 2 :))",), < min(n,m), 

the result is immediate. □ 

Theorem 2.8. If 5J is trace class, L{z; J) can be extended to a contin¬ 
uous function on £\{ —1,1} with 

(2.46) |L( 2 ; J)| < exp |c ||5J||i + ||(5J||i |2 - 1|“^ I 2 + 1|“^ | 

for a universal eonstant, c. //(2.38) holds, L{z;J) ean he extended to all of D 
with 

r r ^ 

(2.47) |L( 2 ;; J)| < expi c 1 + ^ n |a„ - 1| + |6„| 

L n=l 

for a universal constant, c. 

Proof. This follows immediately from (2.22), (2.23), (2.25), and the last 
two lemmas and their proofs. □ 

While we cannot control ||C'^/^(Jo — E{z))~^C^^‘^\\i for arbitrary 2 with 
l^l ^ 1, we can at the crucial points ±1 if we approach along the real axis, 
because of positivity conditions. 

Lemma 2.9. Let C be a positive diagonal trace elass operator. Then 

(2.48) hm (1 - \x\)\\C^/\Jo - E{x))-^C^/^h = 0. 

I^ITl 

X real 
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Theorem 2.10. 


(2.50) 


limsup (1 — |x|) log |L(x; J)| < 0. 

I^ITl 

X real 


Proof. Use (2.30) and (2.18) to write 

L{x; J) = det(l + UC^/‘^{Jo - E{x))-^C^/'^) 
and then (2.15) and (2.31) to obtain 

log|L(x;J)| < \\UC^/\.h-E{x))-^C^/‘^\\i 

< ^\C^/‘^{.h-E{x))-^C^I%. 

The result now follows from the lemma. □ 

Next, we want to find the Taylor coefficients for L{z; J) at z = 0, which 
we will need in the next section. 


(2.51) 


Lemma 2.11. For each fixed h > 0 and \z\ small, 

h 


log 1 - 


Eiz) 




n=l 


where Tn{x) is the Chebyshev polynomial of the first kind: T„(cos6*) = 
cos(n0). In particular, T 2 n+i{ 0 ) = 0 and T 2 n{ 0 ) = (-1)”. 


Proof. Consider the following generating function; 


(2.52) 


g{x, 2 :) = E ^ 


n=l 


The lemma now follows from 
2x 


log 


1 - 


z + z 


-1 


= 2 ( 5 ( 0 , z) - g{x, 2 )] = E i '^nio) - Tn{x) 


by choosing x = hj^. The generation function is well known (Abramowitz and 
Stegun [1, Formula 22.9.8] or Szego [60, Equation 4.7.25]) and easily proved; 
for 0 G M and \z\ < 1, 


dz 


1 


(cos0,z) = - cos{n6)z 

'y • ^ 


n=l 
00 


2z 


Z (.V")” + he-")’ 


n=l 

cos( 0 ) + 2 ; 

2 ^ — 22 cos 6 + 1 

_1 E 

^ dz 


log[l — 2 x 2 + 2 ^] 
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at X = cos6. Integrating this equation from z = 0 proves (2.52) for x G [—1,1] 
and l^l < 1. For more general x one need only consider 9 G C and require 
l^l < exp{ —I Im6<|}. □ 


Lemma 2.12. Let A and B he two self-adjoint m x m matrices. Then 


(2.53) 

where 

(2.54) 


log det {^A — E{z)^ (^B — E{z)^ 


-1 


CX) 


Y,Cn{A,B)z^ 

n=0 


CniA,B) = -^Tr [rn(iA) -T„(is) 


Proof. Let Ai,..., Am be the eigenvalues of A and , hm the eigen¬ 

values of B. Then 


det 


-A-E{z)- 

_B-E{z)_ 


log det 


'A-E{z)' 

B-E{z)_ 



-\j-E{zy 
-Pj - E{z)_ 


log[l - Xj/E{z)] - log[l - Hj/E{z)] 
i=i 


so (2.53)7(2.54) follow from the preceding lemma. 


□ 


Theorem 2.13. IfdJ is trace class, then for each n, T„(J/2) — r„(Jo/2) 
is trace class. Moreover, near z = 0, 


(2.55) 


log[L(z; J)] = ^Cn(J) 2 :” 

n=l 


where 

(2.56) 

In particular, 







(2.57) ci(J) 

(2.58) C2(J) 


-Tr(J-Jo) = -^6m 

m=l 

oo 

-iTr(j2-j2) = -i ^[6^ + 2(a2 

m=l 


!)]• 


Proof. To prove Tn{J/2) — T„( Jo/2) is trace class, we need only show that 
jm_jm _ dJ is trace class, and that’s obvious! Let 5Jn-,F be 

dJn;F extended to £^(Z_|_) by setting it equal to the zero matrix on £‘^{j > n). 
Let Jo,n be Jq with On+i set equal to zero. Then 

^Jn-Akn - EA ^ dJ{Jo - E)-^ 
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in trace norm, which means that 

(2.59) 

K'JQ,n;F ^\^) / 

This convergence is uniform on a small circle about z = 0, so the Taylor series 
coefficients converge. Thus (2.53)/(2.54) imply (2.55)/(2.56). □ 

Next, we look at relations of L{z] J) to certain critical functions beginning 
with the Jost function. As a preliminary, we note (recall is dehned in 
( 2 . 2 )), 

Proposition 2.14. Let SJ be traee class. Then for each z G D\{ — 1, 1}, 

(2.60) lim L(z; J(”)) = 1 

n—^OO 

uniformly on compact subsets of D\{—1,1}. //(2.38) holds, (2.60) holds uni¬ 
formly in z for all z in D. 

Proof. Use (2.16) and (2.24) with B = 0 and the fact that ^ 0 

in the estimates above. □ 

Next, we note what is essentially the expansion of det(J —^^(z)) in minors 
in the first row; 

Proposition 2.15. Let 5J be trace class and z £ l)\{—1,1}. Then 

(2.61) L{z- J) = {E{z) - bi)zL{z- - ajz^L{z- J^^)). 


Proof. Denote by that is, the nx n matrix formed by rows 

and columns A: + 1,..., /c + n of J. Then expanding in minors, 

(2.62) detiE - Jn-,F) = {E- 6i) det(P; - J^^l^.p) - a? det(P; - 

Divide by det(i? — Jo-n-,F) and take n ^ oo using (2.59). (2.61) follows if one 
notes 

det(P/ Jo\n—j\F^ _^ 

det(£' - Jo;n;F) 


by (2.10). 

We now define for z G D\{—1,1} and n = 1,..., oo, 

:”L(z; 


(2.63) 


(ft-.) 

(2.64) 

e 

o 

TP 

II 

(n-i) 


□ 
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Un is called the Jost solution and uq the dost function. The infinite product 
of the a’s converges to a nonzero value since Uj > 0 and J2j juj — 1| < oo. We 
have: 

Theorem 2.16. The Jost solution, Un{z;J), obeys 

(2.65) -\- (pYi E{^Z^']Vjyi ttn^n+1 — 0? ^ — f; • • • 

where ao = 1. Moreover, 

(2.66) lim z~"'Un(z; J) = 1. 

n—^oo 


Proof. (2.61) for J replaced by reads 

L{Z-, J(")) = {E{z) - bn+i)zL{z-, J("+^)) - al+^z^L{z-, j("+2)), 

from which (2.65) follows by multiplying by z^(Jlj^n+i Equation (2.66) 

is just a rewrite of (2.60) because limn^oo O^n Uj = 1- El 

Remarks. 1. If (2.38) holds, one can define Un for z = ±1. 

2. By Wronskian methods, (2.65)/(2.66) uniquely determine Un{z; J). 

Theorem 2.16 lets us improve Theorem 2.5(iii) with an explicit estimate 
on the degree of L{z; J). 

Theorem 2.17. Let 5.J have range n, that is, aj = I if j > n, bj = 0 
if j > u. Then uo(z; J) and so L{z-, J) is a polynomial in z of degree at most 
2n — l. If bn / 0, then L{z-, J) has degree exactly 2n—1. If bn = 0 butun-i / 1, 
then L{z-, J) has degree 2n — 2. 

Proof. The difference equation (2.65) can be rewritten as 


(2.67) 


where 

( 2 . 68 ) 


Un—l 
. Un 


{,E bn)/Un—l Unjan—l j ( Un 

1 ^ J 

1 


Ur 


V^n+l 


An{z) = 


+ 1 — bnZ —ttnZ 
ttn-lZ 0 


If 6J has range n, = Jq and Un = 1- Thus by (2.63), ug{z; J) = z^ if £ > n. 
Therefore by (2.67), 
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( n 

zn+1 

— (cil • • • Ai(^z) ■ ■ ■ An{z') 

= (ai ■ ■ ■ an-i)~^Ai{z) ■ ■ ■ An-liz) 

Since Aj{z) is a quadratic, (2.69) implies uq is a polynomial of degree at most 
2{n — 1) + 1 = 2n — 1. The top left component of Aj contains z^ while 
everything else is of lower order. Proceeding inductively, the top left component 
of Ai{z) ■ ■ ■ An-iiz) is ^2^-2 + 0{z^^-^). Thus if bn / 0, 

uo = -(ai... + 0{z^^~‘^), 


1 - bnZ^ 


O-n—lZ 



,. 0 .) (:j 


proving uq has degree 2n — 1. If 6„ = 0, then 


A„-i(z) 



so inductively, one sees that 




tto = (ai . . . On-l) ^(1 - a 2 _^) 2 ; 2 ” 2 Q('^2n 3^ 


and Uo has degree 2n — 2. 


□ 


Remark. Since the degree of u is the number of its zeros (counting mul¬ 
tiplicities), this can be viewed as a discrete analog of the Regge [46]-Zworski 
[64] resonance counting theorem. 

Recall the definitions (1.2) and (1-15) of the m-function which we will 
denote for now by M{z-, J) = [E{z) — J)il ■ 


Theorem 2.18. 

(2.70) 

(2.71) 


If 5J G Xi then for \z\ < 1 with L{z; J) / 0, we have 


M{z;J) 


zL{z;J^^'>) 
L{z;J) 
ui{z] J) 
uo{z] J) ■ 


Proof. (2.71) follows from (2.70) and (2.63)/(2.64). (2.70) is essentially 
Cramer’s rule. Explicitly, 

M{z-,.J) = lim {E{z) - Jn-,F)n 

n —^oo 

det(E-J«p^) 

= lim ---^ 

n^oo det(£' — Jn-,F) 

= lim WnXmVn 

n^csD 
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where (by (2.59) and (2.10)) 


Wn 




Vn 


detjE - 

det(£^ — Jo;n-l;F) 
det(£^ - Jo;n;F) _ 
det(£^ - Jn-,F) 
det(£^ — Jo;n-l;F) 
det(£' - Jo;„;f) 


--L(z;J«) 


z. 


□ 


Theorem 2.18 allows us to link |tio| and \L\ on \z\ = 1 to Im(M) there; 


Theorem 2.19. Let SJ be trace class. Then for all 6* / 0, vr, the boundary 
value limr|i M{re^^; J) = M(e*®; J) exists. Moreover, 


(2.72) 

Equivalently, 

(2.73) 


no(e*";J) 


2 


Im M(e*^; J) = sin 6. 


L{E^-,J) 


^ Im M(e*^;J) = 



sin0. 


Proof. By (2.64), (2.73) is equivalent to (2.72). If \z\ = I, then E{z) = 
E{z) since z = z~^. Thus, Un{z;J) and Un{z;J) solve the same difference 
equation. Since z~^Un{z; J) ^ 1 and ^ 1, we have that 

(^n[Un{^Z, J'jUn-^-il^Z, T) Un(^Z, J'jUn-^-ii^Z, T)] > Z Z 

Since the Wronskian of two solutions is constant, if z = e*^, 

an[un{e~"^; J)un+i{e'^; J) - Un(e*®; J)ttn+i(e“*^; J)] = 2zsin6'. 

Since oq = 1 and Un{z; J) = Un{z', J), we have that 

(2.74) Im[rto(e*®; J) tti(e*®; J)] =sin0. 

(2.74) implies that wo(e*^; J) / 0 if 6* / 0, tt, so by (2.71), M{z; J) extends 
to T>\{—1,1}. Since ui(e*®; J) = uo{e'^^; J)M{e'^^; J) (by (2.71)), (2.74) is the 
same as (2.72). □ 


If J has no eigenvalues in M\[—2,2] and (2.38) holds so uo{z;J) has a 
continuation to D, then 

/ 1 |■2■K „i9 2 \ 

(2.75) uo{z-,J) = log uo{e'’^] J) d9j 


f2n 


471 Jo — z 


log 


Im M(e*^J) 


sm ( 


de 


(2.76) 


exp 
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(2.77) 

(2.78) 


exp 


1 ^ 
47r Jo e*® — z 


log 


'7r/(2 cos 6 )' 
|sin0| 



{1 - z^) D{z)-^ 


where D is the Szego function defined by (1.29) and f{E) = In the above, 
(2.75) is the Poisson-Jensen formula [48]. It holds because under (2.38), uq is 
bounded on D and by (2.74), and the fact that ui is bounded, log(rio) at worst 
has a logarithmic singularity at ±1. (2.76) follows from (2.72) and (2.77) from 
(1.18). To obtain (2.78) we use 


(1 - z2)2 = exp(^^ I log[sin2 6] doj 


which is the Poisson-Jensen formula for ^(1 


z^)^ if we note that 


(1 _ e-2ie)2 


= 4sin^ 6. 

As a final remark on perturbation theory and dost functions, we note how 
easy they make Szego asymptotics for the polynomials; 


Theorem 2.20. Let J be a Jacobi matrix with SJ trace class. Let Pn{E) 
be an orthonormal polynomial associated to J. Then for \z\ < 1, 

(2.79) lim z'^Pniz -h z"^) = 

n^oo (1 — z^) 

with convergence uniform on compact subsets of D. 


Remarks. 1. By looking at (2.79) near z = 0, one gets results on the 
asymptotics of the leading coefficients of Pn{E), that is, an,n-j in Pn{E) = 
121=0 (^n,kE^; see Szego [60]. 

2. Alternatively, if Qn are the monic polynomials, 

(2.80) lim z'^Qniz + z~^) = 

n^oo (1 — 2;2) 

Proof. This is essentially (2.59). For let 

(2.81) QniE) = det(£' - Jn-p)- 
Expanding in minors in the last rows shows 

(2.82) Qn{E) = {E — bn)Qn-l{E) — a^_iQn-2{E) 

with Qo{E) = 1 and Qi{E) = E — bi. It follows QniE) is the monic orthogonal 
polynomial of degree n (this is well known; see, e.g. [3]). Multiplying (2.81) 
by (oi,... ,a„_i)“^, we see that 

PniE) = (ai... an)~^QniE) 


(2.83) 
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obeys (1.4) and so are the orthonormal polynomials. It follows then from (2.59) 
and (2.9) that 

L{z] J) = lim ^ ^ Qn{z) = lim (1 - z^)z^Qn{z) 

n^oo 2:—W+J-j n—>-oo 

which implies (2.80) and, given (2.83) and lim„^oo(ai • • • Rn)~^ exists, also 
(2.79). □ 


3. The sum rule: First proof 


Following Flaschka [17] and Case [6], [7], the Case sum rules follow from 
the construction of L{z; J), the expansion (2.55) of log[L( 2 ;; J)] at z = 0, the 
formula (2.73) for L(e*^; J) , and the following standard result; 

Proposition 3.1. Let f{z) be analytie in a neighborhood of D, let 
zi,... ,Zm be the zeros of f in D and suppose /(O) / 0. Then 

(3.1) log|/(0)| = — / log /(e*®) de + J2^og\zj\ 

2tt Jo ^ 

and for n = 1,2,..., 


(3.2) Re(a,,) = - r"log|/(e*")| cos(n0) dO — Re 


vr Jo 


m ^—n _ -n 


E 

i=i 


n 


where 

(3.3) 

for \z\ small. 


log 


fiz) 


L/(0)J 


= E 


anZ 


n=l 


Remarks. 1. Of course, (3.1) is Jensen’s formula. (3.2) can be viewed as 
a derivative of the Poisson-Jensen formula, but the proof is so easy we give it. 

2. In our applications, f{z) = f{z) so a„ are real and the zeros are real 
or come in conjugate pairs. Therefore, Re can be dropped from both sides of 
(3.2) and the ” dropped from Zi. 


Proof. Dehne the Blaschke product. 


B{z) 


n 


. , Zj 

J=l J 


Zj-Z 
1 - ZZj 
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While this suffices for the basic sum rule for finite range 5J, which is the 
starting point of our analysis, we note three extensions; 

(1) If f{z) is meromorphic in a neighborhood of D with zeros zi,... ,Zm 
and poles pi,... ^Pk, then (3.1) and (3.2) remain true so long as one makes the 
changes: 


(3.5) 

(3.6) 


i=i 


m —n _ -n 
i 3 


3=1 i=i 

m —n _ k ^ 

\ -^ '^j '^j \ -^ Pj Pj 


-T.- 

i=i 


3 ^3 


for we write f{z) = fi{z)/ Y\j^i{z — Pj) and apply Proposition 3.1 to /i and 
to Y^j=i{z — Pj)- We will use this extension in the next section. 

(2) If / has continuous boundary values on dD, we know its zeros in D 
obey X)^i(l “ kil) < oo (so the Blaschke product converges) and we have 
some control on — log |/(re*®)| as r t 1, one can prove (3.1)-(3.2) by a limiting 
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argument. We could use this to extend the proof of Case’s inequalities to the 
situation '^n\\an — 1| + |&n|] < oo. We first use a Bargmann bound (see [10], 
[19], [20], [27]) to see there are only finitely many zeros for L and (2.73) to 
see the only place log \L\ can be singular is at ±1. The argument in Section 9 
that sup,, /[log_ L{re^^) ]‘^ dO < oo lets us control such potential singularities. 
Since Section 9 will have a proof in the more general case of trace class 6J, 
we do not provide the details. But we would like to emphasize that proving 
the sum rules in generality Case claims in [6], [7] requires overcoming technical 
issues he never addresses. 

(3) The final (one might say ultimate) form of (3.1)/(3.2) applies when 
/ is a Nevanlinna function, that is, / is analytic in D and 

(3.7) sup — / log+ f{re^^) dO < oo, 

0<r<l 271 JQ 

where log_,_(x) = max(log(x), 0). If / is Nevanlinna, then ([48, pg. 311]; essen¬ 
tially one uses (3.1) for f{z/r) with r < 1), 

OO 

(3.8) ^(l-lzjl)<oo 

i=i 

and ([48, pg. 310]) the Blaschke product converges. Moreover (see [48, pp. 247, 
346]), there is a finite real measure on dD so 

(3.9) log f{re^^) dO ^ dij!'^\6) 
weakly, and for Lebesgue a.e. 9, 

(3.10) lim log f{re^^) = log /(e*®) 

rjl 

and 

(3.11) = log /(e*®) d9 + dfiif\9) 

( f) ( f) 

where diJ,s (9) is singular with respect to Lebesgue measure d9 on dD. dfXs (0) 

is called the singular inner component. 

By using (3.1)/(3.2) for f{z/r) with r | 1 and (3.9), we immediately have: 

Theorem 3.2. Let f be a Nevanlinna function on D and let {zj}^^i 
{N = 1,2,..., or oo) be its zeros. Suppose /(O) / 0. Let log /(e*^) be the 
a.e. boundary values of f and d/Us (9) the singular inner eomponent. Then 

1 /*27r 1 r27T 

(3.12) logl/(0)l = ^/o + 
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and for n = 1,2,..., 

1 

(3.13) Re(an) = — log /(e*^) cos{n6) dO 

TT Jo 

1 r ^ z~^ — 

H-/ cos(n0) — Re — - — 

TT Jo n 

where an is given by (3.3). 

We will use this form in Section 9. 

Now suppose that 6J has finite range, and apply Proposition 3.1 to L{z; J). 
Rs zeros in D are exactly the image under E ^ z oi the (simple) eigenvalues 
of J outside [—2,2] (Theorem 2.5(ii)). The expansion of log[L( 2 ;; J)] at z = 0 
is given by Theorem 2.13 and log L(e*®; J) is given by (2.73). We have thus 
proven: 

Theorem 3.3 (Case’s Sum Rules: Finite Rank Case). Suppose 5J has 
finite rank. Then, with |/3i(J)| > \P 2 {J)\ >■■■>! defined so that T 
are the eigenvalues of J outside [—2,2], we have 

1 p2f7V / ' ^ \ 

irl ( Im M(e«») ) -^^ = E ‘°g - E 

^ ^ ^ ^ J n=l 

1 /*27r / sin ^ \ 

( 3 , 15 ) --f Iog(^^-^yos{n0)d0 

= - T”) + I Tt(T„(lj) - T„(lj„)). 

j ^ ^ 


In particular, 
(3.16) 

1 /•27r 


Im M 


-Ode + 


- 7 X! + h X! G'(an), 



= — 1 — log(a^) 

-log\p\% e = 


P + P 1/3] > 1. 


Remarks. 1. Actually, when 5J is finite rank all eigenvalues must lie 
outside [—2, 2] —it is easily checked that the corresponding difference equation 
has no (nonzero) square summable solutions. While eigenvalues may occur at 
—2 or 2 when 5J £ 2i, there are none in (—2,2). This follows from the fact 
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that linirii J) exists for 9 G (0, vr) (see Theorem 2.19) or alternately 

from the fact that one can construct two independent solutions 
whose linear combinations are all non-L^. 


2. In (2.73), 


log|L 


llog 


sin0 
T m M 


OO 

n=l 


the log term is constant and so contributes only to Cq because 

cos(n0) (16 = 0. 

3. As noted, P 2 is Cq + ^C 2 - 

4. We have looked at the combinations of sum rules that give sin^ 6 and 
sin® 9 hoping for another miracle like the one below that for sin^ 9, the function 
G and F that result are positive. But we have not found anything but a mess 
of complicated terms that are not in general positive. 


P 2 is especially useful because of the properties of G and F : 


Proposition 3.4. The function G{a) = — I — 21og(a) for a G (0, 00 ) 

is nonnegative and vanishes only at a = 1. For a — 1 small, 

(3.19) G(a) = 2(a-l)2 + 0((a-l)3). 


Proof. By direct calculations, G(l) = G'{1) = 0 and 


G"{a) = 


2(1+ a2) 


> 2 


so G{a) > (a —1)^ (since G(l) = G'(l) = 0). (3.19) follows from G"(l) = 4. □ 


Proposition 3.5. The function F{e) given by (3.18) is positive through¬ 
out its domain, {|e| > 2}. It is even, inereases with inereasing |e|, and for 
|e| — 2 small, 

(3.20) F{e) = |(|e| - 2)^/2 + 0((|e| - 2)^). 

In addition, 

(3.21) F{e) <l{e^ 


Proof. Let R{/3) = — /? 2 _ log |/3| ) for /3 > 1 and compute 

This shows that R{f3) is increasing. It also follows that 


R'{P) = 2iP-2f + 0{iP-lf) 
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and since /3>1, (/3 + l)//3<2 and (3 ^ < 1 so 

<2(/3-1)2. 

As i?(l) = 0, we have 

(3.22) R{f3) <§(/?- 
and 

(3.23) i?(/3) = |(/3_l)3 + o((^_l)4). 

Because F{—e) = -F(e), which is simple to check, we can suppose e > 2 so 
(3 > 1. As/3 = i[e + — 4 ] is an increasing function of e, F{e) = -R(/3) is 

an increasing function of e > 2. Moreover, (3 — 1 = {e — 2)^/^ + 0(e — 2) and 
so (3.23) implies (3.20). Lastly, 

{(3-l)<(3-f3-^ = ^e‘^-A, 

so (3.22) implies (3.21). □ 


4. The sum rule: Second proof 


In this section, we will provide a second proof of the sum rules that never 
mentions a perturbation determinant or a Jost function explicitly. We do 
this not only because it is nice to have another proof, but because this proof 
works in a situation where we a priori know the m-function is analytic in a 
neighborhood of D and the other proof does not apply. And this is a situation 
we will meet in proving Theorem 6. On the other hand, while we could prove 
Theorems 1, 2, 3, 5, 6 without Jost functions, we definitely need them in our 
proof in Section 9 of the Co-sum rule for the trace class case. 

The second proof of the sum rules is based on the continued fraction 
expansion of m (1.5). Explicitly, we need, 

(4.1) -M(z; J)-^ = -{z + z-^) + 6i + alM{z; 


which one obtains either from the Weyl solution method of looking at M (see 
[24], [56]) or by writing M as a limit of ratio of determinants 


det(E(z) - 


(4.2) 


M{z-, J) = lim 


det{E{z) - Jn-p) 


and expanding the denominator in minors in the first row. For any J, (4.1) 
holds for z € D. Suppose that we know M has a meromorphic continuation to 
a neighborhood of D and consider (4.1) with z = e*®: 

(4.3) -M(e*^ J)-^ = -2cos0 + 6i + a?M(e*®; jW). 
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Taking imaginary parts of both sides, 

oid. 


( 4 . 4 ) 




or, letting 


{Note: because 
(4.5) M{z] J) = {z + z~^ — J)ii = z{l + z^ 

near zero, g is analytic in D), we have 

sin0 / sin0 

To see where this is heading. 



g{z;J) = 


M{z-J) 


zJ)il = z + 0(z^} 


= logai+log g{e^^-J) 


Theorem 4.1. Suppose M{z;J) is meromorphic in a neighborhood of D. 
Then J and have finitely many eigenvalues outside [—2, 2] and if 

1 /•27r / SlU f) \ ^ 

{with f3j as in Theorem 3.3), then 

(4.8) C'o(J) = -log(ai) + Co(J«). 

In particular, if 5J is finite rank, then the Co sum rule holds: 


(4.9) C'o(J) = - ^ log(a„). 

n=l 


Proof. The eigenvalues, Ej, of J outside [—2,2] are precisely the poles of 
m{E; J) and so the poles of M{z; J) under E^ = z^ + zj^. By (4.1), the poles 
of M{z; J^^)) are exactly the zeros of M{z; J). Thus {/3j(J)“^} are the poles of 
M{z;J) and are its zeros. Since g{0;J) = 1 by (4.5), (3.1)/(3.5) 

becomes 


1 


log(5(e*^J) 


de 


^log(|/3,(J)|)+^log(/3,(J«) 


(4.6) and this formula imply (4.8). By (4.3), if M{z;J) is meromorphic in 
a neighborhood of D, so is M{z; So we can iterate (4.8). The free M 

function is 


( 4 . 10 ) 


M{z-, Jo) = z 
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(e.g., by (2.7) with m = n = 1), so C'o(Jo) = 0 and thus, if 5J is finite rank, 
the remainder is zero after finitely many steps. □ 


To get the higher-order sum rules, we need to compute the power series 
for log( 5 '(z; J)) about z = 0. For low-order, we can do this by hand. Indeed, 
by (4.1) and (4.5) for 

g{z;J) = {z[{z + z~^) - bi - afz + 0{z‘^)])~^ 

= {1 - hz - (al - l)z‘^ + Oiz^))-^ 

= l + biz + i{ai-l) + bl)z^ + 0{z^) 
so since log(l + w) = w — + 0{w^), 

(4.11) log{g{z; J)) = biz + {^bl + aj- l)z‘^ + O(z^). 

Therefore, by mimicking the proof of Theorem 4.1, but using (3.2)/(3.6) in 
place of (3.1)/(3.5), we have 


Theorem 4.2. Suppose M{z;J) is meromorphic in a neighborhood of D. 
Let 


(4.12) 

Cn{J) 



sin0 \ 
ImM(e*^) / 


cos(n0) -|- 






Then 


(4.13) Ci(J) = 6i+Ci(jd)) 

(4.14) C 2 (J) = [\hl + {al-l) +C 2 (/^))]. 

If 

1 ^277 / Sin 0 \ 

(4.15) = + 

with F given by (3.18), then writing G{a) = of — 1 — 21og(a) as in (3.17) 

(4.16) P2{J) = I bl +I G(ai) + P2{J^^^). 

In particular, ifSJ is finite rank, we have the sum rules Ci, C 2 , P 2 o/(3.15)/(3.16). 


To go to order larger than two, we expand \og{g{z; J)) systematically as 
follows: We begin by noting that by (4.2) (Cramer’s rule). 


(4.17) 
where 

(4.18) 


g{z;J) = lim gn{z;J) 

n —>00 


9n{z;J) 


z ^ det( 2 ; + z ^ 

det(2: -|- z~^ — Ju-f) 

1 det(l-j^(z)-Vl^2i;,.) 

1-1-5:^ det(l — £'(z)“^ Jn;E) 


(4.19) 
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where we used z{E{z)) = l + z'^ and the fact that because the numerator has a 
matrix of order one less than the denominator, we get an extra factor of E{z). 
We now use Lemma 2.11, writing Fj{x) for j[rj(0) — Tj{x/2)], 

(4.20) 

\og gn{z]J) = 


OO r 

log(l + ^2) + E Tr {e,{/X,f)) - Tr (i",(Jn;F)) 


j=i 


(4.21) 




= -log(l + z2)-^ —(-1) 


-1J 




j=l 3 


where we picked up the first sum because Jn-F has dimension one greater than 
J^li-F so the Tj{0) terms in Fj{Jn-p) and contribute differently. Notice 

^ — (-1)^ = -log(l + z2) 

j=l 3 

so the first two terms cancel! Since gn{z] J) converges to g{z; J) in a neighbor¬ 
hood of z = 0, its Taylor coefficients converge. Thus 


Proposition 4.3. For each j. 


(4.22) 


a 


:(J, J*'^^)= lim 






exists, and for z small, 


(4.23) 


logg{z; J) = ^ — aj{J, 
j=i 7 


Remark. Since 




m+1 


if m > 7, the difference of traces on the right side of (4.22) is constant for 
n > j, so one need not take the limit. 

Plugging this into the machine that gives Theorem 4.1 and Theorem 4.2, 
we obtain 


Theorem 4.4. Suppose M{z;J) is meromorphic in a neighborhood of D. 
Let Cn{J) he given by (4.12) and a by (4.22). Then 

(4.24) Cn{J) = - an{J, J(^)) + Cn{J^^^). 

n 

In particular, if 5J is finite rank, we have the sum rule Cn o/(3.15). 
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Proof. The only remaining point is why if 5J is finite rank, we have re¬ 
covered the same sum rule as in (3.15). Iterating (4.24) when J has rank m 
gives 

ly m 

(4.25) Cn{J) = 

= Tr Tn{\jt,F) - Tn{\Jo,f.-m.-,F) 

while (3.15) reads 

(4.26) Cn{J) = ^Tr r„(ij) -T„(iJo) 

2 ' 

= Tr - Tr . 

That (4.25) and (4.26) are the same is a consequence of Proposition 2.2. □ 

5. Entropy and lower semicontinnity 
of the Szegd and quasi-Szegd terms 

In the sum rules Cq and P 2 of most interest to us, there appear two terms 
involving integrals of logarithms; 

(5-1) Z{J) = ^ r^logf n ) 

^ ^ ^ ’ All Jo ^VImM(e*®, J)/ 

and 

One should think of M as related to the original spectral measure on 
ct(J) D [—2, 2] as 

(5.3) ImM(e^^) = tt (2cos^) 


in which case, (5.1), (5.2) can be rewritten 



Our main result in this section is to view Z and Q as functions of n and to 
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prove if ^ weakly, then Z{^n) (resp. Q{Hn)) obeys 

(5.6) Z{p) < liminf Q{l^) < hminf Q(^n)) 

that is, that Z and Q are weakly lower semicontinuous. This will let us prove 
sum rule-type inequalities in great generality. 

The basic idea of the proof will be to write variational principles for Z and 
Q as suprema of weakly continuous functions. Indeed, as Totik has pointed 
out to US, Szego’s theorem (as extended to the general, not only a.c., case [2], 
[18]) gives what is essentially Z{J) by a variational principle; explicitly, 

(5.7) exp{i/^''log(!^)<i«} = mf[i£ <He)\ 

where P runs through all polynomials with P{0) = 1, which can be used 
to prove the semicontinuity we need for Z. It is an interesting question of 
what is the relation between (5.7) and the variational principle (5.16) below. 
It also would be interesting to know if there is an analog of (5.7) to prove 
semicontinuity of Q. 

We will deduce the semicontinuity by providing a variational principle. 
We originally found the variational principle based on the theory of Legendre 
transforms, then realized that the result was reminiscent of the inverse Gibbs 
variation principle for entropy (see [55, pg. 271] for historical remarks; the 
principle was first written down by Lanford-Robinson [35]) and then realized 
that the quantities of interest to us aren’t merely reminiscent of entropy, they 
are exactly relative entropies where /x is the second variable rather than the 
first one that is usually varied. We have located the upper semicontinuity of 
the relative entropy in the second variable in the literature (see, e.g., [12], [34], 
[44]), but not in the generality we need it, so especially since the proof is easy, 
we provide it below. We use the notation 

(5.8) logj_(x) = max(± log(x), 0). 


Definition. Let ^, 1 / be finite Borel measures on a compact Hausdorff 
space, X. We define the entropy of ^ relative to u, S{iJ, \ n), by 


(5.9) 


S{n I v) 


—oo if is not u-ac 


Remarks. 1. Since log_(x) = log_,_(x ^) < x ^ and 

I{^) # o}) < KM < oo 

the integral in (5.9) can only diverge to —oo, not to -|-oo. 
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2. If d/i = / du, then 

(5.10) S{^\u) = - J /log(/)du, 
the more usual formula for entropy. 

Lemma 5.1. Let ^ he a probability measure. Then 

(5.11) S{ijL I v) < logv'{X). 

In particular, if v is also a probability measure, 

(5.12) S{ix I u) < 0. 

Equality holds in (5.12) if and only if fi = n. 

Proof. If fj, is not u-ac, (5.11)/(5.12) is trivial, so suppose p, = f du and 
let 


(5-13) dh = X{x\f(x)^Q}dv 

so V and p. are mutually ac. Then, 

(5.14) S{p, I v) = 




< log 


(5.15) 


ydnJ 
dv 
dfi 

= logu(X) 

< logu{X) 


dp, 


where we used Jensen’s inequality for the concave function log(x). For equality 
to hold in (5.12), we need equality in (5.15) (which says v = v) and in (5.14), 
which says, since log is strictly convex, that du/dp is a constant. When v{X) = 
p{X) = 1, this says v = p. □ 


Theorem 5.2. For all p, u. 


(5.16) 


S{p \ v) = inf 


F{x) du 


+ logF) dp{x) 


where the inf is taken over all real-valued continuous functions F with min^, GxFi. 

> 0 . 


Proof. Let us use the notation 

G{F,p,n)= J F{x) dn - J {I + log F) dp{x) 
for any nonnegative function F with F G L^{dv) and logT G L^{dp). 
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Suppose first that ^ is p-ac with dfi = f dv and F is positive and contin¬ 
uous. Let A = {x \ f{x) / 0} and define D by (5.13). As log(a) is concave, 
log(a) < a — 1 so for a, 6 > 0, 

(5.17) ab~^ > 1 + log(a6“^) = 1 + log(a) — log6. 

Thus for X G A, 

F{x)f{x)~^ > 1 + logF(x) — log/(x). 

Integrating with dfj, and using 

J F{x)dv> J F{x)dv = J F{x)f{x)~^dix, 

we have that 

J F{x) du > y (1 + log F(x)) d^{x) + S{fi\ u) 
or 

(5.18) S{fi \ u) < G{F, ^1,1/). 

To get equality in (5.16), take F = f so J dfi and J F du cancel. Of course, 
/ may not be continuous or strictly positive, so we need an approximation 
argument. Given N,£, let 


r N if f{x) > N 
fN,e{x) = < /(x) if e < /(x) < N 
[ e if /(x) < e. 

Let fe,N,eix) be continuous functions with s < fe,N,e < N so that as £ ^ oo, 
fe,N,e fN,e in L^(A, dfi + du). For IV > 1, < 1 -t- /, so we have 

- J ^Og{fN,e) dfi = ^(^) / 

< „(Miog[i + hfi] <oo 

and thus, since — log(/Ar^£) increases as e J, 0, fN,e=o = linieio fN,e has 
log/7v,£=o £ L^{dfi) and the integrals converge. It follows that as £ ^ oo 
and then e | 0, 


^ ^ Q{fN,e=0, 

We now take N ^ oo. By monotonicity, — J log fN,e=o dfu. converges to 
— / log / djj, which may be infinite. In addition, f fN,e=o du — J dfj, ^ 0 so 
^(/a,£= 0 )/n, p) ^ S{fj, I r), and we have proven (5.18). 
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Next, suppose /i is not u-ac. Thus, there is a Borel subset A C X with 
fi(A) > 0 and v{A) = 0. By regularity of measures, we can find K <Z A 
compact and for any e, Ue open so K C A C Ue and 

(5.19) > 0 u(t4) < e. 

By Urysohn’s lemma, find continuous with 

(5.20) 1 < F^{x) < £~^ all X, F^ = £~^ on K, F^ = 1 on X\Ue- 
Then 

j F,du< uiX\Ue) + £-^l^{Ue) < V{X) + 1 

while 

J{1+log Fs)d^ > log{e~^)^{K) 
so 

Q{Fe,ii, v) < v{X) + 1 - \og{£~^)^{K) -oo 

as e I 0, proving the right side of (5.16) is —oo. □ 

As an infimum of continuous functions is upper semicontinuous, we have 

Corollary 5.3. S{^ \ u) is jointly weakly upper semieontinuous in p, 

and u, that is, if Pn —^ P o,nd v, then 

S{p\v)> lim sup S'(/Un I Vn). 

n 

Remarks. 1. In our applications, pn will be fixed. 

2. This proof can handle functions other than log. If J log{{dn/dp)~^) dp 
is replaced by / G{{du/dp))dp where G is an arbitrary increasing concave 
function with limj^^o G{y) = oo, there is a variational principle where 1 + logT 
in (5.18) is replaced by H{F{x)) with H{y) = m.lx{xy — G{x)). 


To apply this to Z and Q, we note 

Proposition 5.4. 

(a) Let 


(5.21) 

dpQ{E) = ^ ^4 - dE. 

ZTT 

Then 


(5.22) 

Q{J) = -\S{pQ 1 pj). 

(b) Let 


(5.23) 

, , 1 dE 

dpi(E) =- . 

^ ^ ^ vr V4-E2 
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Then 

(5.24) Z{J) = -i log(2) - i 5(^1 I MJ). 


Remarks. 1. Both hq and fii are probability measures, as is easily checked 
by setting E = 2 cos 6. 

2. dfiQ is the spectral measure for Jq. For M{z]Jo) = z and thus 
ImM(e*®; Jq) = sin0 so m{E; .Jq) = ^\/A — E"^ and ^ ImmdE = d^Q. 

3. dfii is the spectral measure for the whole-line free Jacobi matrix and 
also for the half-line matrix with = 0, oi = \/2, 02 = 03 = • • • = 1. An 
easy way to see this is to note that after E = 2cos0, dni{9) = ^d9 and 
so the orthogonal polynomials are precisely the normalized scaled Chebyshev 
polynomials of the first kind that have the given values of aj. 

Proof, (a) Follows immediately from (5.16) if we note that 
duo _ ^ 

dfx dE / dE 2 tt dfis^/dE 


(b) As above, 


Thus 

where 

c 


^ = 2(4 - E^)-^ _ 

dfi 2 tt d^a.c/dE 


Z{J) ^ c - I iij), 


- I 

2-K J-2 

r 27 r 


log 


4-.F2 


\/4 - F;2 dE 


1 

— / log[2sin^ 0] 

471 Jo 

1 1 

-log(2)-F— / log|sin6»| (i6l 
2 271 Jo 

1 1 f27r I _ Jd 

-log{2) + -f log^-d« 

1 log(2) + log(l) = -1 log(2) 


where we used Jensen’s formula for f{z) = ^(1 ~ z"^) to do the integral. 


□ 


Remark. As a check on our arithmetic, consider the Jacobi matrix J with 
ai = \/2 and all other a’s and 6’s the same as for Jq so d^j is d^i. The sum 
rule, Co, for this case says that 

Z(J) = -log(v/2) = -llog2 
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since there are no eigenvalues and ai = \f2. But ni = /xj, so 5(/ii | /xj) 
= 0. This shows once again that c = — ^log2 (actually, it is essentially the 
calculation we did—done the long way around!). 

Given this proposition, Lemma 5.1, and Corollary 5.3, we have 

Theorem 5.5. For any Jacobi matrix, 

(5.25) Q{J) > 0 
and 

(5.26) Z{J) > -i log(2). 

U FJn f'J weakly, then 

(5.27) Z(J) < liminf Z(J„). 
and 

(5.28) Q{J) < liminf Q{Jn)- 

We will call (5.27) and (5.28) lower semicontinuity of Z and Q. 

6. Fun and games with eigenvalues 

Recall that Jn denotes the Jacobi matrix with truncated perturbation, as 
given by (2.1). In trying to get sum rules, we will approximate J by Jn and 
need to estimate eigenvalues of Jn in terms of eigenvalues of J. Throughout this 
section, X denotes a continuous function on M with X{x) = X{—x), X(x) = 0 
if |x| < 2, and X is monotone increasing in [2,oo). Our goal is to prove: 

Theorem 6.1. For any J and all n, we have N^{Jn) < N^{J) + 1 and 

(i) \Et{Jn)\<\Et{J)\ + l, 

(ii) |F±^i(Jn)| < |F±(J)|. 

In particular, for any function X of the type described above, 

N±{Jn) Ar±(J) 

(6.1) ^ XiEf{Jn))<XiEt{J) + l)+ 

j=i j=i 

Theorem 6.2. If J — Jq is compact, then 

N±(E) N±(j) 

(6-2) E X{Ef{Jn))= E 

j=i i=i 

This quantity may be infinite. 

Proof of Theorem 6.1. To prove these results, we pass from J to Jn in 
several intermediate steps. 
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(1) We pass from J to Jn-F- 

(2) We pass from to Jn-F i dn,n = ^n-F '''^here dn^n is the matrix with 1 
in the n, n place and zero elsewhere. 

(3) We take a direct sum of J^.p and Jq ± di^i. 

(4) We pass from this direct sum to Jn- 

Step 1. Jn-F is just a restriction of J (to ^^({1, • • • ,?^}))- The min-max 
principle [45] implies that under restrictions, the most positive eigenvalues 
become less positive and the most negative, less negative. It follows that 

(6.3) N^{Jn-,F) < N^{J) 

(6.4) ±EfiJn-,F) < ±Ef{J). 

Step 2 . To study E^, we add dnp, and to study Ej , we subtract dn,n- 
The added operator dn,n has two critical properties; It is rank one and its 
norm is one. From the norm condition, we see 

(6.5) \Et{J^^p)-EtiJn-,F)\<l 
so 

(6.6 ) ^i+(4f) ^ Ei{Jn-,F) + l 

< El{J) + l. 

(Note (6.5) and ( 6 . 6 ) hold for all indices j, not just j = 1, but we only need 
j = 1.) Because dn,n is rank 1, and positive, we have 

^m+li'^n-,F) < E^pi{J^.p) < E^{Jn-^F) 

and so, by (6.4), 

(6.7) E+^,iJ+p)<E+{J) 
and thus also 

( 6 . 8 ) iV±(J±^)<iV±(J) + l. 

Step 3. Take the direct sum of J^.p and Jq ± dn. This should be inter¬ 
preted as a matrix with entries 

f k,i<n 

^n;F ® ('^0 ® ^ll) ^ ^ \ {Jq ® n,£—n ^ 71 

[ 0 otherwise. 

Since Jq ± dn has no eigenvalues, (6.7) and ( 6 . 8 ) still hold. 
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Step 4. Go from the direct sum to Jn- In the + case, we add the 2x2 
matrix in sites n,n + 1: 

^ -1 n 


and, in the — case. 


dJ+ = 


dJ- = 


1 -1 


1 1 
1 1 


dJ~^ is negative, so it moves eigenvalues down, while dJ is positive. Thus 

< E++i{J+f) < E+{J) 

and 




□ 


Proof of Theorem 6.2. We have, since J — Jq is compact, that 

> 0 . 


Thus 

(6.9) 


\Jn — <7|| < sup \bm\ + 2 sup \a„ 

m>n+l m>n 




If ^{Ef{J)) = oo, then, by (6.9), for all fixed m, 

N^(Jn) m 

lim inf ^ XiEfiJn)) > liminf ^X(£;±(J„)) 
j=i i=i 

m 

= ^x(E±(j)) 
i=i 

so taking m to infinity, (6.2) results. 

If the sum is hnite, (6.9), dominated convergence and (6.1) imply (6.2). □ 


7. Jacobi data dominate spectral data in P 2 

Our goal in this section is to prove Theorem 5. Explicitly, for a Jacobi 
matrix, J, let 

CXD OO 

02(j) = J !:'>? +5 Eooj) 


(7.1) 
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with G = a? — 1 — 21og(a) as in (3.17). For a probability measure, ^ on M, 
define 

(7.2) ft{„) = i «<*« + £ PiEj) 

where Ej are the mass points of jj, outside [—2,2] and F is given by (3.18). 
Recall that Imil£j(e*^) = irdfis^c/dE at E = 2cos0. We will let /Uj be the 
measure associated with J by the spectral theorem and the Jacobi matrix 
associated to fj,. 

Then Theorem 5 says 

Theorem 7.1. If J — Jq is Hilbert-Schmidt so that D 2 {J) < oo, then 
P2{hj) < oo and 

(7.3) P2{iaj) < D2{J). 

Proof. Let Jn be a truncation of J given by (2.1). Then D2{Jn) is mono¬ 
tone increasing with limit D2{J). This is finite because J — Jq is Hilbert- 
Schmidt. By the definition (5.5), 

(7.4) P2{fr) = Q(J)+^F(E£ 

j 

= -^^S{^^Q,^i)+J2F{E,) 

j 

by (5.22). Since Q > 0 (5.25) and F > 0 (Proposition 3.5), (7.4) is a sum 
of positive terms. Moreover, by Theorem 6.2, J2j F{Ej{.Jn)) J2j P{PjiJ)) 
even if the right side is infinite. As Jn ^ J in Hilbert-Schmidt sense, {Jn—E)~^ 
converges (in norm) to (J — E)~^ for all F G C \ M. This implies that fij^ 
converges weakly to /a and so by (5.28), Q{J) < limsup Q{Jn)- R follows that 

F2 (Jm) < limsup g(Jn) + ^F(F£J„)) 

= limsup Il 2 (J„) (by Theorem 3.3) 

= D2{J). 

Thus P 2 {jj) < oo and (7.3) holds. □ 

The result in this section is essentially a quantitative version of the main 
result in Deift-Killip [13]. 

8 . Spectral data dominate Jacobi data in P2 

Our goal in this section is to prove the following, which is essentially 
Theorem 6: 
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Theorem 8.1. If fj, is a probability measure with P2{h) < oo, then 

( 8 . 1 ) D2{J^) < P2iti) 

and so is Hilbert-Schmidt. 

The idea of the proof is to start with a case where we have the sum rule and 
then pass to successively more general cases where we can prove an inequality 
of the form (8.1). There will be three steps; 

(1) Prove the inequality in the case is meromorphic in a neighborhood 
oiD. 

(2) Prove the inequality in the case p, > dfiQ where 5 is a positive real number 
and fiQ is the free Jacobi measure (5.18). 

(3) Prove the inequality in the case ^ 2 ( 11 ) < co. 

Proposition 8.2. Let J be a Jacobi matrix for which has a mero¬ 
morphic continuation to a neighborhood of D. Then 

( 8 . 2 ) D2{J) < P2{J). 

Proof. By Theorem 4.2, 

P2{J) = Ibj + ^G{ai) + P2{J^^^). 

so iterating, 

m m 

P2iJ) = lJ2b] + hT.Giaj) + P2{J^^^) 

i=i i=i 

m m 

j=l j=0 

since P 2 (J(™)) > 0. Now G > 0, so we can take m ^ 00 and obtain (8.2). □ 

Remark. If has a meromorphic continuation into {z \ \z\ < rj} for 
some rj > 1, then by a theorem of Geronimo [21], \ an — 1\ p"' + \bn \ p"' <00 

for all p < p, so the sum rule also follows from the methods of Section 3. We 
prefer to avoid the use of Geronimo’s theorem. 

Given any J and associated M-function M{z;J), there is a natural ap¬ 
proximating family of M-functions meromorphic in a neighborhood of D. 

Lemma 8.3. Let he the M-function of a probability measure p obeying 
the Blumenthal- Weyl condition, and define 

(8.3) M^P{z) = r~^M^{rz) 

for 0 < r < 1. Then, there is a set of probability measures so that = 
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Proof. Return to the E variable. Since (z) is meromorphic in a neigh¬ 
borhood of D with > 0 if Imz > 0, 

m^^\E) = -M^''\z{E)) 

(where z{E) -|- z{E)~^ = E with \z\ < 1) is meromorphic on C\[—2,2] and 
Herglotz. It follows that it is the Borel transform of a measure of total 
weight liuiE^oo—Em^'^^E) = liuizio z~^ Mjf'\z) = 1. □ 

Proposition 8.4. Let /i be a probability measure obeying the Blumen- 
thal- Weyl condition and 

(8.4) p. > (5^0 

where /uq is the free Jaeobi measure {the measure with Mfj,g{z) = z) and 5 > 0. 
Then 

(8.5) D2{Jff) < P2(^)- 


Proof. We claim that 


(8.6) limsup / —log ImM (r)(e*^) d6< —log ImM^(e*^) 


rfl 


de. 


Accepting (8.6) for the moment, let us complete the proof. The eigenvalues of 
that lie outside [—2,2] correspond to /3’s of the form 




h{J) 


for those k with |/3fc(J)| < r. Thus F{E^{J^(r))) is monotone increasing to 

E{E^{Jf^)), so (8.6) shows that 


(8.7) 


P 2 {h) > limsup P 2 {ijS'^^). 


Moreover, M^(r){z) M^{z) uniformly on compact subsets of D which 

means that the continued fraction parameters for vn^'^ (E), which are the Jacobi 
coefficients, must converge. Thus for any N, 


N N-l 


i X! + 7 X! 

i=i i=i 


< 


< 

< 


N 


N-l 


lim 

rfl 


i=i 


{r)\ 


J=1 

lim inf D 2 {J^(t) ) 

lim inf P 2 {pS^'^) (by Proposition 8.2) 
P 2 {h) (by (8.7)) 


so (8.5) follows by taking N 


oo. 
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Thus, we need only prove ( 8 . 6 ). Since M^(r){9) = r ^M^(re*®) ^ M^(e*®) 
for a.e. 6, Fatou’s lemma implies that 

( 8 . 8 ) liminf J log+ ImM^(r)(e*®) dO > J log+|ImM^( 6 *)| dO. 

On the other hand, (8.4) implies |ImM^( 2 ;)| > ( 5 |Im 2 ;|, so ImM^(r)( 2 :) > 
5 |Im 2 :|. Thus uniformly in r, 

(8.9) ImM^(j.) (e*^) >(5|sin0|. 

Thus 

log_ ImM^(r)(e*®) < — log(5 — log |sin 0 | , 
so, by the dominated convergence theorem, 

lim J log_( ImM^M ) d9 = J log_(|ImM^( 0 )|) 

This, together with ( 8 . 8 ) and — log(x) = — log^(x) + log_(x) implies ( 8 . 6 ). □ 

Remark. Semicontinuity of the entropy and (8.9) actually imply one has 
equality for the limit in ( 8 . 6 ) rather than inequality for the lim sup. 

Proof of Theorem 8.1. For each 5 £ (0,1), let = (1 — 5)fi + Since 
Us obeys (8.4) and the Blumenthal-Weyl criterion, 

( 8 . 10 ) D2{J^g) < P2{h&)- 

Let Ms = and note that 

Im Ms{e^^) = (1 — d) Im M(e*^) + 5 sin 9 
so 

log \mMs{e^^) = log(l — d) + log ImM(e*®) + --- sind 

1 — 0 

We see that up to the convergent log(l —d) factor, log Im.Ms{e'^^) is monotone 
in (5, so by the monotone convergence theorem, 

(8.11) P2(/J) = lim F)2(/i5) 

(the eigenvalue terms are constant in 8, since the point masses of ^s have the 
same positions as those of ^!). 

On the other hand, since ^s ^ T weakly, as in the last proof. 


(8.12) D 2 {Jff) < liminf D 2 (J/,J. 

(8.10)-(8.12) imply (8.1). 


□ 
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9. Consequences of the Cq sum rule 

In this section, we will study the Cq sum rule and, in particular, we will 
prove Nevai’s conjecture (Theorem 2) and several results showing that control 
of the eigenvalues can have strong consequences for J and //j, specifically 
Theorems T and 7. While Nevai’s conjecture will be easy, the more complex 
results will involve some machinery, so we provide this overview: 

(1) By employing semicontinuity of the Szego term, we easily get a Co-in¬ 
equality that implies Theorems 2 and 7 and the part of Theorem A' that 
says J — Jo is Hilbert-Schmidt. 

(2) We prove Theorem 4.1 under great generality when there are no eigen¬ 
values and use that to prove a semicontinuity in the other direction, and 
thereby show that the Szego condition implies a Co-equality when there 
are no eigenvalues, including conditional convergence of J2ni^n — !)• 

(3) We use the existence of a Co-equality to prove a Ci-equality, and thereby 
conditional convergence oiJ2n^n- 

(4) Returning to the trace class case, we prove that the perturbation deter¬ 
minant is a Nevanlinna function with no singular inner part, and thereby 
prove a sum rule in the Nevai conjecture situation. 

Theorem 9.1 (= Theorem 3). Let J be a Jacobi matrix with cJess(J) C 
[—2, 2] and 

(9.1) 

(9.2) 

Then 

(i) CTess(J) = [- 2 , 2 ]. 

(ii) The Szego condition holds; that is, 

Z{J) < oo 

with Z given by (5.1). 

(hi) crac(J) = [— 2 , 2 ]; indeed, the essential support of a^c is [— 2 , 2 ]. 


^efc(J)^/^ < oo, 
k 

N 

lim sup log(aj) > —oo. 
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Remarks. 1. We emphasize (9.2) says > —oo, not < oo, that is, it is a 
condition which prevents the a^s from being too small (on average). 

2. We will see below that (9.1) and (9.2) also imply \aj — 1| ^ 0 and 
\bj\ ^0 and that at least inequality holds for the Cq sum rule: 

N 

(9.3) Z{.J) < ^log|/3fc(J)| - limsup^log(aj) 

k ^ j=i 

holds. 


Proof. Pick A^i, • • • (tending to oo) so that 


(9.4) 


inf 

e 


Ne 

Illog(ai) 

j=i 


> —oo 


and let be given by (2.1). By Theorem 3.3, 

Ne 

(9.5) Z{Jm^) < -^log(aj)+ ^log(|/3fc(JArJ|) 
i=i 

Ne 

- + + 2^os(l/3i(^)| + 2) 

i=i 


where in (9.5) we used Theorem 6.1 and the fact that the /3 solving ei(J) +1 = 
/3 + (i.e., 1 + (5i-\- = (i + (i~^) has (3 < /3i(J) + 2. For later purposes, 

we note that if \hn{J)\ + \an{J) — 1| ^ 0, Theorem 6.2 implies we can drop the 
last term in the limit. 

Now use (5.27) and (9.5) to see that 


Z{J) < liminf Z{Jn() < oo. 

This proves (ii). But (ii) implies > 0 a.e. on i? G [—2,2], that is, 
[—2,2] is the essential support of ;Uac- That proves (iii). (i) is then immedi¬ 
ate. □ 


Proof of Theorem 2 {Nevafs conjeeture). We need only check that J — Jq 
trace class implies (9.1) and (9.2). The finiteness of (9.1) follows from a bound 
of Hundertmark-Simon [27], 

^[ |e,(J)| \ek{J) + 4| ]i/2 < ^ |6J + 2 |a„ - 1| 

n 

where ek{J) = |F'^| — 2 so |e| |e + 4| = — 4. 

Condition (9.2) is immediate for, as is well-known, aj > 0 and 

X)(|aj| —1) < oo implies is absolutely convergent, that is, |log(fli)l < oo. 

□ 
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Corollary 9.2 (= Theorem 7). A discrete half-line Schrodinger opera¬ 
tor (i.e., Un = 1) with aess{J) C [—2,2] and < oo has a^c = [—2,2]. 

This is, of course, a special case of Theorem 9.1 but a striking one discussed 
further in Section 10. In particular, if = 1 and bn = n~‘^Wn where a < ^ and 
Wn are identically distributed independent random variables with distribution 
g{X) dX with g G and supp( 5 ) bounded, then it is known that [—2, 2] is 
dense pure point spectrum (see Simon [54]). It follows that J must also have 
infinitely many eigenvalues outside [—2,2], indeed, enough that 
= oo. 

Next, we deduce some additional aspects of Theorem 4h 

Corollary 9.3.//(Tess(«4) C [—2,2] and (9.1), (9.2) hold, then J—JqG Z 2 , 
that is, 

(9.6) ^6^+^(an - 1)2 < 00 . 

Proof. By Theorem 6, (9.6) holds A < 00 , and Q{J) (given 

by (5.21)) is finite. By (9.1) and ek{J)^^^ < ei{J)ek{jy^‘^ , we have that 
X]efc(J)^/2 < 00 . Moreover, Z{J) < 00 (i.e.. Theorem 9.1) implies Q{J) < 00 . 
For, in any event, flmM dO < 00 implies 



Q{J) < 00 . 


What remains to be shown of Theorem 4' is the existence of the conditional 
sums. We will start with X](on — !)• Because Y^(an — 1)^ < 00 , it is easy 
to see that — 1) is conditionally convergent if and only if X]log(an) is 

conditionally convergent. By (9.5) and the fact that J — Jq is compact, we 
have: 


Proposition 9.4. If (9.2) holds and cr{J) C [—2,2], that is, no eigen¬ 
values outside [—2,2], then 


(9.7) 


Z{J)< 


lim sup 


N 

J2^og{aj) 


'-1=1 
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We are heading towards a proof that 


(9.8) 


Z{J)> 


lim inf 


N 

Illog(aj) 


'-i=i 


from which it follows that the limit exists and equals Z{J). 


Lemma 9.5. If a{J) C [—2,2], then log[ 2 ;“^M( 2 ;; J)] lies in every H'P{D) 
space for p < oo. In particular, J) is a Nevanlinna function with no 

singular inner part. 


Proof. In D\(—1,0), we can define KigM{z-,J) C (—7r,7r) and Argz C 
(—71, tt) since ImM( 2 ;; J)/Imz > 0. Thus g{z;J) = z~^M{z;J) in the same 
region has argument in (—7r,7r). But Arg^i is single-valued and continuous 
across (—1,0) since M has no poles and precisely one zero at z = 0. Thus 
P^Tgg G L“. It follows by Riesz’s theorem on conjugate functions ([48, pg. 351]) 
that log( 5 () G HP{D) for any p < oo. Since it lies in g is Nevanlinna. Since 
for p > 1, any function, F, has boundary values F{re^^) T(e*^) in 

log( 5 ') has no singular part in its boundary value. □ 

Proposition 9.6. Let a{J) C [—2,2]. Suppose Z{J) < oo. Let Co,Cn 
be given by (4.10) and (4.15) {where the j3{J) terms are absent). Then the 
step-by-step sum rules, (4.8), (4.13), (4.14), (4.24) hold. In particular, 

(9.9) Z{J) = -log(ai) + Z(J«) 

(9.10) CiiJ) = 6i + Ci(jW). 


Proof. (4.4) and therefore (4.1) hold. Thus, we only need apply Theo¬ 
rem 3.2 to g, noting that we have just proven that g has no singular inner 
part. □ 

Theorem 9.7. If J is such that Z{J) < oo and (j{J) C [—2,2], then 

(i) limAT^oo Ejli log(aj) exists. 

(ii) The limit in (i) is —Z{J). 

(hi) 

(9.11) lim Z(j(")) =0 (= Z{Jo)) 

n—^oo 

Proof. By (9.9), 

n 

Z(J)+^log(a,) = Z(jW). 
i=i 


(9.12) 
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Since J — Jq G £ 2 , /UJq weakly, and so, by (5.27), liminf > 0, 

or by (9.12), 


(9.13) 


r n 


lim inf 


Illog(aj) 

w=i 


> -Z{J). 


But (9.7) says 


r n 


lim sup 


Illog(ai) 

w'=i 


< -Z{J). 


Thus the limit exists and equals Z{J), proving (i) and (ii). Moreover, by (9.12), 
(i) and (ii) imply (hi). □ 


If Z(-) had a positive integrand, (9.11) would immediately imply that 
Ci(jW) ^ 0 as n ^ 00 , in which case, iterating (9.10) would imply that 
J2]=ibj is conditionally convergent. Z{-) does not have a positive integrand 
but a theme is that concavity often lets us treat it as if it does. Our goal is to 
use (9.11) and the related lim„^oo = 0 (which follows from Theorem 5) 

to still prove that C'i(J’^"'^) ^ 0. We begin with 


Lemma 9.8. 
J fndfj, < 1, and 

(9.14) 

Then 

(9.15) 


Let dn be a probability measure and suppose fn > 0, 


J log{fn) dp, = 0. 

J |log(/n)| dp + J \fn -1\ dp^ 0. 


Proof. Let 

(9.16) H{y) =-log{y)-l + y. 

Then 


(i) H{y) >0 for all y. 

(ii) inf|y_i|>ei7(|/) > 0. 

(hi) H{y) > iy if y > 8. 

(i) is concavity of log(y), (ii) is strict concavity, and (hi) holds because 
— logy — 1 + ^y is monotone on (2, 00 ) and > 0 at y = 8 since log(8) is slightly 
more than 2. 

Since J (fn — ^)dp < 0, (9.14) and (i) implies that 

J fn(x) dp(x) 1 


(9.17) 
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(9.18) 


lim [ H{fn{x))dfj,{x) 0. 

n—^oo J 


Since H >0, (ii) and the above imply ^ 1 in measure: 
(9.19) fi{{x I \fn{x) - 1| > e}) ^ 0. 

By (i), (hi) and (9.18), 


(9.20) / 

Jfn{x)>8 

Now (9.19)/(9.20) imply that 


\fn{x)\ dfi 0 . 


J \fn{x) - 1| dfi{x) 0 


and this together with (9.18) implies / |log(/„)| dfi = 0. 

Proposition 9.9. Suppose Z{J) < oo and (t{J) C [—2,2]. Then 


(9.21) 


lim / log 


ImM(e*®, 


de = 0 . 


Proof. By (9.11), the result is true if | • | is dropped. Thus it suffices to 


lim / log_ 


d9 = 0 


n-lloy_/-^-VlmM(e*^ JH); " 

or equivalently, 

(9.22) lim = 

V Sm0 / 

Now, let dpQ{6) = ^sin^ 9 d9 and fn{9) = (sin0)“^ ImM(e*®, By 

(1.19), 


d9 = 0. 


(9.23) 


fni9)dfio{9) < 1 


and by Theorem 5 (and Corollary 9.3, which implies || — Jolli ^ 0)) 

J ^og{fn{9)) dfio{9) 0 

so, by Lemma 9.8, we control |log| and so log^; that is, 

.. r,.... = 


(9.24) lim [ log, f- 


Thus, to prove (9.22), we need only prove 

(9.25) lim lim sup/ , (In. 


d9 = 0. 
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To do this, use 

< log+(a)+log_(6) = 21og+(a^/^)+ log_(6) 

< 2a^/2 ^log_(6) 

with a = sin0ImM(e*^, and b = sin^ 0. The contribution of log_(6) in 
(9.25) is integrable and n-independent, and so goes to zero as e J, 0. The 
contribution of the 2a^/^ term is, by the Schwartz inequality, bounded by 

1 /2 

(4e)V2(^4 f_jn{e)di,^{e)^ 

also goes to zero as e J, 0. Thus (9.25) is proven. □ 

The following concludes the proofs of Theorems 4 and 4'. 


Theorem 9.10. If Z{J) < oo and a{J) C [—2,2], then 


(9.26) 


N 

lim 22 bj exists and equals 
N^oo ^i 
J=1 


1 \ 
271 7o ^ VimTf (e*®)/ 


cos(0) dO. 


Proof. By Proposition 9.9, C'i(J*^"')) ^ 0 and, by (9.10), 

n 

Ci{J) = J2bj + CiiJ^^2- □ 

i=i 


As a final topic in this section, we return to the general trace class case 
where we want to prove that the Cq (and other) sum rules hold; that is, we 
want to improve the inequality (9.5) to an equality. The key will be to show 
that in this case, the perturbation determinant is a Nevanlinna function with 
vanishing inner singular part. 


Proposition 9.11. Let J — Jq be trace class. Then, the perturbation 
determinant L{z] J) is in Nevanlinna class. 

Proof. By (2.19), if is given by (2.1), then 

(9.27) L{z-Jn)^L{z-J) 

uniformly on compact subsets of D. Thus 


(9.28) 

sup 

0<r<l 



log+ L(re*'^; 


J) 


de 

— < 
277 “ 


^27r 


sup sup / log+ L(re ;J„) 

n 0<r<lJ0 ' 

/•27r I 

JO 


sup / log+ L(e* ;J„) 

n JO 


de 

2^ 


de 

2^ 
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where (9.28) follows from the monotonicity of the integral in r (see [48, pg. 336]) 
and the fact that L(z; Jn) is a polynomial. 

In (9.28), write log+ |L| = logjLj + log_ |L|. By Jensen’s formula, (3.1), 
and L(0; J) = 1, 

f2ir JQ 

/ log L(e*®;Jn) 73-=-X^log|/ 3 j(Jn)I 

Jo 271 ^ 

and this is uniformly bounded in n by the ^ Lieb-Thirring inequality of 
Hundertmark-Simon [27], together with Theorem 6.1. On the other hand, 
by (2.73), 

(9.29) 21„g_ |L(e»; J„)| = log_ ( R 4 

n—1 

< 2 ^ log_(aj) + 21og(sin0) + log.,. (^ImM(e*^; Jn) sin0^ 
i=i 

since log_ |a6/c| = [log(a) + log(6) — log(c)]_ < log_(a) + log_(6) + log^.(c). 
The first term in (9.29) is 0-independent and uniformly bounded in n since 
[uj — 1| < oo. The second term is integrable. For the final term, we note 
that log^.(y) < y so by (1.19), the integral over 9 is uniformly bounded. □ 

Remark. Our proof that L is Nevanlinna used X]fc(^fc( J))^^^ < oo as input. 
If we could find a proof that did not use this a priori, we would have, as a 
consequence, a new proof that 6fc( J)^'^^ < oo since X][l “ /Jfc( J)~^] < oo is 
a general property of Nevanlinna functions. 

Proposition 9.12. If5J GXi, the singular inner part of L{z;J), if any, 
is a positive point mass at z = 1 and/or at z = —1. 

Proof. By Theorem 2.8, L{z] J) is continuous on D\{—1, 1} and by (2.73), 
it is nonvanishing on {e*^ | 0 7 ^ 0,7r}. It follows that on any closed interval, 
I C (0,7r) U (7r,27r), log L{re^^,J) dO converges to an absolutely continuous 
measure, so the support of the singular inner part is {±1}. 

Returning to (9.29) and using log^.(x) = 2 log^.(a;^/^) < , we see that 

log_ |L| lies in that is, 

sup sup /[log_ , J„) ]^ ^ < oo 

n 0<r<lJ 271 

and this implies log_ L{re^^,J) has an a.c. measure as its boundary value. 
Thus ±1 can only be positive pure points. □ 
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Remark. We will shortly prove L has no singular inner part. However, we 
can ask a closely related question. If |X]log(a„)| < oo and < oo 

so Z{ J) < oo, does the sum rule always hold or is there potentially a positive 
singular part in some suitable object? 

Theorem 2.10 will be the key to proving that T(x; J) has no pure point 
singular part. The issue is whether the Blaschke product can mask the polar 
singularity, since, if not, (2.46) says there is no polar singularity in L which 
combines the singular inner part, outer factor, and Blaschke product. Experts 
that we have consulted tell us that the idea that Blaschke products cannot mask 
poles goes back to Littlewood and is known to experts, although our approach 
in the next lemma seems to be a new and interesting way of discussing this; 

Lemma 9.13. Let f{z) be a Nevanlinna function on D. Then for any 
00 e dD, 

(9.30) lim [log(l - r)"^]“^ [ log /(ye*^“) dy = 2/r^({6'o}). 

fTi Jo 
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so 

B{z) =Y[b^.{z) 

Zi 

and note that for 0 < x < 1 and any a E -D, 

1 > \ba{x)\ > b\c,\{x) , 
so 

(9.33) 0<-log|B(x)| <y]-log b\^.\{x) . 

Zi 

Thus, also without loss, we can suppose all the zeros Zi lie on (0,1). 

If a G (0,1), a straightforward calculation (or Maple!) shows 

(9.34) [ -log| 6 „(x)| dx = alog("-) + ^ logf-^—V 

JO \a/ a \1 — a/ 

We claim that for a universal constant C and r > |, a > ^, 

(9.35) — / log | 6 q,(x)| dx < (7(1 — a) log(l — r)“^. 

Jo 

Accepting (9.35) for the moment, by (9.33), we have for r > |, 



where r]{a) is the right side of (9.34). Dividing by log(l — r) ^ and using 
T]{a) < oo, we see 

limsupf-—log|S(x)| dxj] < (7 (1 “ «i)- 

Llog(l-r) U JJ 

Taking n ^ oo, we see that the limsup is 0. Since — log \B{x)\ > 0, the limit 
is 0 as required by (9.32). Thus, the proof is reduced to establishing (9.35). 

Note first that if 1 > a > ^, ^(1 + a) < 4. Moreover, if g{a) = alog(^), 
then g"{a) = — ^ < 0, so g{a) < (1 — a). Thus, if r/(a) is the right side of 

(9.34) and a > ^, then 

(9.36) 17 ( 0 ) < 1 — a + 4(1 — a) log 



Suppose now 

(9.37) 


1 — a > (1 — r)^. 


1 

1 — a 
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Then 


rr 

(9.38) — / log|6a(x)| dx < ri{a) < (1 — a) 

Jo 


< (1 - a) 


1 + 8 log 


1 + 4 log 
1 


1 — a 


1 — r 


by (9.37). 

On the other hand, suppose 
(9.39) (l-a)<(l-rf. 

By an elementary estimate (see, e.g., [48, pg. 310]), 


(9.40) 


|1 - br,(x)l < 


1 — X 


(1 - a). 


If (9.39) holds and x < r, then, by (9.40), 

2(1 - r )2 


(9.41) 


|1 - 6n,(a:)| < 


1-. S2(l-r)<- 


since r is supposed larger than |. IfttG(^,l), then 

— logu= [ — < 2(1 — u) 
Ju y 


so if (9.41) holds. 


- log(6a(x)) < 2(1 - ba{x)) < 


4(1 — a) 
1 — x 


and so 

(9.42) [ — log(6a(x)) dx < 4(1 — a) log(l — r)“^. 

do 

We have thus proven (9.38) if (9.37) holds and (9.42) if (9.39) holds. Together 
this proves (9.35). □ 

Theorem 9.14. Let J — Jq be trace class. Then the Nevanlinna func¬ 
tion L{z] J) has a vanishing singular inner component and all the sum rules 
Co, Cl,... hold with no singular term. 

Proof. By Proposition 9.12, the only possible singular parts are positive 
points at ±1. By (9.30) and the estimate (2.50), these point masses are absent. 
Thus the singular part vanishes and the sum rules hold by Theorem 3.2. □ 
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10. Whole-line Schrodinger operators with no bound states 

Our goal in this section is to prove Theorem 8 that the only whole-line 
Schrodinger operator with a{W) C [—2, 2] is Wq, the free operator. We do this 
here because it illustrates two themes: that absence of bound states is a strong 
assertion and that sum rules can be very powerful tools. 

Given two sequences of real numbers {&n}^-oo with a„ > 0, 

we will denote by W the operator on defined by 

(10.1) iWu)n = an-lUn-l + bnUn + ttn+lUn+l- 

Wo is the operator with On = 1, 6n = 0. The result we will prove is; 

Theorem 10.1. Let W be a whole-line operator with an ^ 1 and a{W) C 
[—2, 2]. Then W = Wo, that is, bn = 0. 

The proof works if 

r m 1 

limsup ^ log(aj) > 0. 

m—*00 j — — ^ 

The strategy of the proof will be to establish analogs of the Cq and C 2 sum 
rules. Unlike the half-line case, the integrand inside the Szegd-like integral will 
be nonnegative. The Co sum rule will then imply this integrand is zero and 
the C 2 sum will therefore yield J2n = 0. As a preliminary, we note: 

Proposition 10.2. Ifon^l and a{W) C [—2,2], then 
particular, bn ^ 0 as \n\ -^ 00 . 

Proof. Let J be a Jacobi matrix obtained by restricting to {1, 2,...}. By 
the min-max principle [45], (t(J) C [—2,2]. By Corollary 9.3, J2'?^=ibn < 00 . 
Similarly, by restricting to {0, —1, —2,...,}, we obtain X)n=-oo < 00 . □ 

Let for n = 1,2,... be the operator with 

[ f ’" ‘ 

(10.2) bp = bj if \j\ < n 

[ bP = 0 if |j| > n. 

Then, Proposition 10.2 and the proofs of Theorems 6.1 and 6.2 immediately 
imply: 

Theorem 10.3. If an ^ I and a{W) C [—2,2], then W^'^^ has at most 
four eigenvalues in M\[—2,2] {up to two in each of {— 00 , —2) and (2, 00 )) and 
for j = 1,...,4, 

(10.3) |ej(lT(”))| ^ 0 


as n 


00 . 
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Note. As in the Jacobi case, ej{W) is a relabeling of E^{W) — 2 in 
decreasing order. 

To get the sum rules, we need to study whole-line perturbation determi¬ 
nants. We will use the same notation as for the half-line, allowing the context 
to distinguish the two cases. So, let 5W = W — Wq be trace class and define 

( 10 . 4 ) L{z;W) = det{{W - E{z){Wo - E{z))-^) 

( 10 . 5 ) = dei{l+ 5W{Wo - E{z))-^) 

where as usual, E{z) = z + z~^. 

The calculation of the perturbation series for L is algebraic and so imme¬ 
diately extends to imply: 

Proposition 10.4. If6W is trace class, for each n, Tn{W/2) — Tn{Wo/2) 
is trace class. Moreover, for \z\ small, 

OO 

( 10 . 6 ) log[L{z-,W)] = '£cn{W)z^ 

n=l 

where Cn{W) is 

( 10 . 7 ) Cn{W) = Tv(Tn[^w) -Tn[^Wo)y 
In particular, 

OO 

( 10 . 8 ) c2{W) = --J2bl + 2{al-l). 

m=l 

The free resolvent, (ITo — E{z))~^, has matrix elements that we can com¬ 
pute as we did to get (2.9), 

(H'o - = -(2-> - zj-'zl”-”! 

which has poles at z = ±1. We immediately get 

Proposition 10.5. If 5W is finite rank, L{z] J) is a rational function 
on C with possible singularities only at z = ±1. 

Remarks. 1. If 5W has bo = 1, all other elements zero, then 

Hz., w) = i- (z-' - z)-' = 

has poles at ±1, so poles can occur. 

2. The rank one operator 


R{z) = -(^-1 
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is such that if SW = with C finite rank, then 

C^/\Wo - E{z)-^ - R{z)]C^/^ 

is entire. Using this, one can see L(z; J) has a pole of order at most 1 when 
5W is finite rank. We will see this below in another way. 

If 2 : G l)\{—1,1}, we can define a dost solution Un{z;W) so that (2.65) 
holds for all n G Z and 

(10.9) lim z"”u+( 2 ;IU) = 1. 

n —^00 

Moreover, if 5W has finite rank, is a polynomial in 2 ; for each n > 0. 
Moreover, for n < 0, z~'^u^ is a polynomial in 2 by using (2.67). 

Similarly, we can construct u~ solving (2.65) for all n G Z with 

lim z^u-{z;W) = 1. 

n^—oo 

As above, if SW is finite rank, u~ is a polynomial in 2 if n < 0 and for n > 0, 
z^u~ is a polynomial. 

Proposition 10.6. Let 5W be trace class. Then for z G l)\{ —1,1} and 
all n G Z, 

/ 00 

(10.10) L{z;W) = ( 2 “^- 2 )"M Y[ aj 

^ j=-oo 

an[u+(2;IU)u“+i(2;IU) - u"(2; IU)u++i(2; IP)]. 

Proof. Both sides of (10.10) are continuous in W, so we need only prove 
the result when 5W is finite rank. Moreover, by constancy of the Wronskian, 
the right side of (10.10) is independent of n so we need only prove (10.10) when 
I 2 I < 1 and n is very negative-so negative it is to the left of the support of 
6W, that is, choose ii so On = 1, 6n = 0 if n < —R, and we will prove that 

(10.10) holds for n < —R. 

For n < —R, 2 "' and 2 “” are two solutions of (2.65) so in that region we 

have 

( 10 . 11 ) = aez"^ + /?£ 2 “”. 

Taking the Wronskian of given by (10.11) and u~ = 2 “"' at some point 
n < —R, we see 

/ 00 

(10.12) RHS of (10.10) = aU Y[ aj 

^j=-oo 

Let Wn, Wo-n be the Jacobi matrices on + 1, n + 2,..., }) obtained by 

truncation. On the one hand, as with the proof of (2.74), for \z\ < 1, 

(10.13) L(z;W)= lim det{{Wn - E{z))iWn-o - E{z))-^) 

n—5 -—00 
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and on the other hand, for n < —R, by (2.64), 

RHS of (10.13) = (Y[an\-^u+{z]W). 

^ n ' 

Thus 

(10.14) L{z;W) = lim z"”u+(z;lT) 

\ / n—>—cxD 

\ r) ■' 


m an ] lim {ae + Pe 

' n—*—csD 


—2n\ 


^ n 

since \z\ < 1 and n —oo. Comparing (10.12) and (10.14) yields (10.10). □ 


Note. In (10.11), ae,Pe use ‘T’ for “left” since they are related to scatter¬ 
ing from the left. 

Corollary 10.7. If 6W is finite rank, then (1 — z‘^)L{z; W) is a poly¬ 
nomial and, in particular, L{z-, W) is a rational function. 

Proof. By (10.10), this is equivalent to z{uQuf — u^uf) being a polyno¬ 
mial. But Uq , zuf, Uq , and are all polynomials. □ 

Let 6W be finite rank. Since L is meromorphic in a neighborhood of D 
and analytic in D, Proposition 3.1 immediately implies the following sum rule: 


Theorem 10.8. If 6W is finite rank, then 

N{W) 


(10.15) Co: 

(10.16) Cn: 


for n > 1. 


1 72 ^ 


TZTT '' ^ 

/ log L{P^-,W) de= Y. logl/?.WI 

•^0 


1 Ctt 

vr Jo 


log L(e*®;lT) cos{n6) dO 
N{W) 


^ T 1/5” - an - f - r„(iiF„) 


The final element of our proof is an inequality for L(e*®; W) that depends 
on what a physicist would call conservation of probability. 


Proposition 10.9. Let 5W he trace class. Then for all 6 0 ,tt, 


(10.17) 


L(e*^lT) 


OO 

> n 

j=-oo 
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Proof. As above, we can suppose that bW is finite range. Choose R so 
that all nonzero matrix elements of bW have indices lying within (—i?, i?). By 
(10.12), (10.17) is equivalent to 

(10.18) la^l > 1 

where ai is given by (10.11). 

Since ufiz; W) is real for z real, we have 

ut{P,W) = ut{z-W). 

Thus for z = , 0 / 0, TT, and n < —R, 

u+(e*^lT) = af(e*®)e*’^® +/3f(e*®)e-*^® 

u+{e-^^;W) = . 

Computing the Wronskian of the left-hand sides for n > R, where = z'^ 
and then the Wronskian of the right-hand sides for n < —R, we find 

i(sin6») = f(sin6')[|a£|^ - \(3(\^] 

or, since 0 / 0, vr, 

(10.19) \ai\^ = l + m^ 

from which (10.18) is obvious. □ 

Remark. In terms of the transmission and reflection coefficients of scat¬ 
tering theory [61], = 1/t, = r/t, (10.19) is |rp -|- = 1 and (10.18) is 

\t\ < 1. 

We are now ready for 

Proof of Theorem 10.1. Let be given by (10.2). Then, by (10.3) and 
Co (10.15), 

(10.20) y log L(e*®;IT(”)) = 0. 

Since a„ = 1, (10.17) implies log L(e*®; > 0, and so (10.20) implies 

(10.21) lim — / cos(26')log L(e*®;IT(’")) (10 = 0. 

n^oo 271 Jo 

By (10.8), Om = 1, C 2 , and (10.3), we see 

lim = 0 


which implies 6 = 0. 


□ 
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Finally, a remark on why this result holds that could provide a second 
proof (without sum rules) if one worked out some messy details. Here is a part 
of the idea; 


Proposition 10.10. Let {bn} be a bounded sequenee and W the associ¬ 
ated whole-line Sehrodinger operator {with On = 1). Let 


(10.22) 

yl(a)=^6„e-“W 

If 


(10.23) 

limsup A{a) > 0, 

«i0 

W has spectrum in (2, oo) 

, and if 

(10.24) 

liminf A{a) < 0, 

aiO 

then W has spectrum in ( 

-00,-2). 


Proof. Let (10.23) hold and set ipa{n) = e "1^1/2. Then 


{Woipa){n) 

It follows that 


2cosh(^) if n / 0 

[2cosh(^) — 2sinh(^)](/3Q,(n) if n = 0. 


(10.25) 

(ifa, Wipa){n) 


2 cosh 

2\\ipA 


a 
'2 

+ 2 


+ A{a) - 2sinh^^^ 


cosh(|)-l 


+ A{a) - 2sinh^^y 


Now, sinh(Q;/2) ^ 0 as a J, 0 and since = 0{a~^) and cosh(Q;/2) — 1 = 

O(a^), 2[cosh(a/2) — 1 ]||(/?q,|P ^ 0 as a | 0. If there is a sequence with 
limH(Q;ji) > 0, for n large, {(pa^,iW — 2)ipa^) > 0 which implies there is 
spectrum in (2,oo). 

If (10.24) holds, use (fain) = (—and a similar calculation to 
deduce {ipa„,iW + 2)ipaJ <0. □ 


This proof is essentially a variant of the weak coupling theory of Simon 
[50]. Those ideas immediately show that if 

(10.26) n |6n| < oo 

and J2bn = 0 (so Proposition(lO.lO) does not apply), then W has eigenvalues 
in both (2,oo) and (—oo,2) unless 6 = 0. This reproves Theorem 10.1 when 
(10.25) holds by providing explicit eigenvalues outside [—2, 2]. It is likely using 
these ideas as extended in [4], [31], one can provide an alternate proof of 
Theorem 10.1. In any event, the result is illuminated. 
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Added Notes. During the refereeing process, several results were obtained 
which relate to this paper. In connection with Theorem 10.1, D. Damanik, 
D. Hundertmark, R. Killip, and B. Simon (to appear) have proved that if the 
essential spectrum of a whole- or half-line Schrodinger operator is contained in 
[—2,2], then it is a compact perturbation of the free operator. B. Simon and 
A. Zlatos (to appear) have studied when the Cq sum rule holds, have simplified 
the proofs of Theorems 7.1 and 9.14, and have extended Theorem A' to the case 
where one assumes (1-12) rather than that there is no discrete spectrum. 
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